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and/or  fire  size  hold).  However,  two  special  preexisting  atmospheric  conditions  could  alter 
the  nature  of  the  large-area  fire  significantly,  and  these  are  also  examined  by  simple 
approximate  analyses.  First,  preexistence  of  a  well-developed  atmospheric  cyclone  over  the 
fire  area  could  bring  conservation-of-angular-momentum  considerations  to  paramount 
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importance;  criteria  for  onset  of  a  combustlon-heat-intenslfled  cyclone 
(firestorm)  are  proposed.  Second,  preexistence  of  a  strong  sustained  wind 
of  fairly  well-defined  direction  could  permit  the  fire  to  propagate  rapidly 
beyond  the  area  encompassing  the  Initial  Ignitions,  with  only  an  exhaustion 
of  fuel,  noteworthy  change  of  weather  conditions,  and/or  significant  topo- 
graplc  fire  break  terminating  the  flame  sperad. 
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I.  INTRODUCTION 

The  primary  consequences  of  thermonuclear  weapons  outside  the  area  of 
150-mb  peak  overpressure  are  likely  to  be  of  an  incendiary  nature.  Yet 
examination  of  basic  references  {e.g.,  Glasstone  and  Dolan  1977)  readily 
reveals  the  relatively  miniscule  effort  devoted  to,  and  rather  debatable 
statements  about,  fire  effects;  blast,  radiative,  and  fallout  effects  have 
been  examined  more  extensively  and  are  far  better  quantified. 

It  is  now  recognized  that  a  much  better  understanding  of  urban-scale 
fires  is  needed  if  one  is  to  anticipate  the  effects  on  the  atmosphere  of 
such  events  (Crutzen  and  Birks  1982;  Pyne  1982;  Tilling  1982;  Holden  1983; 
Raloff  1983;  Wooley  and  Bishop  1983).  In  particular,  information  on  plume 
height,  smoke  content  and  early  scavenging  [particularly  fires  burning  in  a 
humid  environment  (Morton  1957)]  deserves  attention. 

In  this  report  we  present  some  preliminary  studies  of  urban-scale 
uncontrolled  fires.  In  these  broad-brush  preliminary  examinations,  the  aim 
is  to  include,  in  the  simplest  reasonable  manner,  only  the  most  important 


2. 


CATEGORIZING  URBAN-SCALE  FIRES 


Incendiary  effects  of  thermonuclear  weapons  on  urban  and  wildland  envi¬ 
ronments  can  involve  multiple  simultaneous  ignitions  in  the  presence  of  heavy 
locating  of  combustibles,  blasted  into  strewn  debris.*  It  seems  plausible 
that  the  urban  environment  after  blast  is  likely  to  have  a  far  more  continuous 
(even  if  still  discrete)  distribution  of  combustibles,  such  that  the  ground- 
level/canopy-level  fuel  distribution  of  forests  is  more  closely  approached  in 
the  blasted  city.  However,  whereas  ladder  (intermediate-height)  fuels  are 
limited  in  forests  to  down  branches,  young  trees,  mixed  species,  lichen,  etc., 
ladder  fuels  should  be  plentiful  in  a  multistory  urban  environment  even  after 
blast  (in  those  outer  areas  where  the  incendiary  aftermath  constitutes  the  prin¬ 
cipal  nuclear  effects).  Subsequent  merger  of  initially  separated  ignitions 
across  normal ly-effective  firebreaks  may  be  abetted  in  the  blasted  environment, 
especially  since  the  disruption  of  transportation,  communication,  and  water- 
distribution,  casualties  in  trained  personnel  and  in  firefighting  apparati ,  etc., 
probably  curtail  effective  countermeasures . 

Attempts  to  anticipate  the  precise  number  and  nature  of  ignitions  in  a  spe¬ 
cific  locale  seems  difficult  in  view  of  the  accidental  character  of  firestarts 
after  earthquakes.  Rather,  the  occurrence  of  ignitions  and  early  spread  is 
accepted  here.  It  seems  more  productive  to  seek  insight  concerning  (1)  the  rate 
and  areal  extent  of  the  spread  of  the  fire  front  through  thinner  (and  more  rapidly 
burned)  fuels;  (2)  the  completeness  of  consumption  cf  thicker  (and  more  slowly 
burned)  fuels;  (3)  the  magnitude  and  overall  character  of  winds;  (4)  the  properties 


Despite  intense  burning  in  the  aftermath  of  the  atomic  bombings  of  Hiroshima  and 
Nagasaki  in  August  1945,  the  suggestion  arose  in  tne  late  1970's  that  blasted- 
gas  arrival  at  the  site  of  radiation-precursor-initiated  diffusion  flames  would 
effect  forced-convection  extinction  of  the  primary  fires.  However,  blowing  out 
of  the  diffusion  flame  supported  by  a  sublimating  or  evaporating  fuel  results 
in  rapid  lowering  of  the  temperature  of  liquid  or  solid  fuel  by  convective  cooling, 
such  that,  even  after  the  wind  decreases,  burning  is  not  restored  because  the 
fuel  vapor  pressure  is  too  low;  thus,  even  a  rather  brief  positive  phase  suf¬ 
fices  to  extinguish  burning  permanently.  However,  wildland  and  urban  combusti¬ 
bles  are  predominantly  char-forming  polymers  that  retain  heat  in  a  porous- 
carbonaceous-matrix  layer  that  forms  under  thermal  degradation  to  envelop  the 
pristine  core;  the  heat- retaining  properties  of  the  char  sustains  pyrolyzation 
of  further  fuel  for  a  finite  interval  of  convective  cooling,  such  that  autoreig¬ 
nition  of  hot  outgassed  fuel  vapor  may  occur  at  the  end  of  the  positive  phase. 

In  fact,  heterogeneous  oxidation  of  the  char  by  the  wind  during  the  positive 
phase  helps  counter  convective  cooling  (Carrier  et  al .  1982). 
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of  the  buoyant  convection  (its  height  and  width,  its  propensity  to  loft  soot 
and  firebrands,  its  tendency  to  Initiate  precipitation);  and  (5)  the  mechanisms 
and  time  scale  for  ultimate  extinction  of  the  fire. 

If  the  primary  (radiation-precursor-initiated)  and  secondary  (blast-caused) 
fires  tend  to  persist  and  merge,  then  one  has  what  is  termed  henceforth  a  big 
(i.e.,  large-area)  fire.  In  that  many  phenomena  may  be  of  significance,  it 
seems  particularly  difficult  to  give  generally  valid  results  via  simple  analyses. 
That  is,  even  if  one  postulates  that  (1)  the  ambient  meteorological  state  (atmos¬ 
pheric  stability  and  moisture  content  and  cloudiness,  and  wind  magnitude  and 
direction,  as  a  function  of  altitude)  is  known;  (2)  the  spatial  distribution  of 
combustibles  (the  total  mass,  the  size  distribution,  ignitabi 1 ity ,  outgassing 
characteristics,  exothermic ity,  char-forming  and  soot- forming  characteristics,  and 
moisture  content)  is  known;  and  (3)  the  to.-ography  is  known,  still  it  is  quite 
challenging  to  describe  the  evolution  of  the  fire  from  a  given  initial  condition. 
However,  work  submitted  for  publication  on  the  necking  in  and  height  of  the  buoy¬ 
ant  plume,  on  the  magnitude  and  mechanism  for  enhanced  entrainment  into. the  lower 
flaming  portion  of  the  plume,  and  on  the  treatment  of  situations  involving  appre¬ 
ciable  buoyancy  persisting  high  into  the  tropopause  is  presented  as  Appendix  A. 

Two  special  atmospheric  circumstances  result  in  a  small  number  of  phenomena 
being  of  particular  importance  in  the  mass,  momentum,  and  heat  balances,  and  thus 
would  seem  to  offer  simple  analysis  better  prospects  for  yielding  insight.  First, 
in  a  relatively  small  number  of  cases,  there  may  happen  to  exist  a  preexisting, 
relatively  weak  cyclone  at  the  site  of  the  large  fire;  as  a  result,  spin-up  under 
convectively- induced  advection  may  already  be  under  way,  and  the  fire  may  lead 
to  a  so-called  firestorm.  Criteria  for  the  onset  of  a  nonpropagating  but  '■'xcep- 
tionally  intense  fire  (in  which  conservation  of  angular  momentum  furnishes  key 
constraints)  have  been  developed  in  other  work  also  submitted  for  publication, 
and  are  presented  here  as  Appendix  B.  Second,  in  a  larger  number  of  cases,  there 
may  exist  a  significant  preexisting  low-altitude  wind,  persisting  in  direction 
and  sustained  in  magnitude.  The  result  is  a  wind-aided  flame  spread  across  strewn 
debris.  Notes  on  this  subject  have  not  been  submitted  for  publication  by  the 
authors  and  constitute  the  remainder  of  the  main  body  of  this  final  report. 


3.  WIND-AIDED  FLAME  SPREAD  ACROSS  STREWN-DEBRIS-TYPE  FUEL 


3. 1  Preliminary  Comments 

Hit  spectacular  urban  fires  of  modern  times  have  usually  been  associated 
with  strong  sustained  winds  (London,  1666;  Lisbon,  1775;  Moscow,  1812;  Chicago, 
1871;  Boston,  1872;  Baltimore,  1904;  Tokyo/ Yokohama,  1923;  Bandon,  Oregon,  1936; 
Tokyo,  1945).  Many  of  the  memorable  wildland  fires  (Pyne  1982)  were  conse¬ 
quences  of  wind-aided  flame  spread  (Miramichi  River  Valley,  New  Brunswick, 

Canada,  1825;  Peshtigo,  Wisconsin,  1871;  Hinkley,  Minnesota,  1894;  northern 
Rockies,  1910;  Cloquet,  Minnesota,  1918;  Matilija,  southern  California,  1932; 
Tillamook  County,  Oregon,  1933;  Shoshone  National  Forest,  Wyoming,  1937;  Vic¬ 
toria,  Australia,  1939;  southern  Maine,  1947;  Inja  and  Santa  Monica  Mountains, 
southern  California,  1956;  Sundance  Mountain,  Idaho,  1967;  Mack  Lake,  Michigan, 
1980;  Victoria,  Australia,  1983).  This  list  is  hardly  exhaustive.  What  it  sug¬ 
gests  is  that  ignition  often  occurs  in  heavily  fuel-laden  areas  in  times  of 
drought,  but  it  is  the  coincidence  of  persistent  winds  of  appreciable  speed  that 
causes  a  “blow-up."  The  arising  of  strong  winds  initiates  a  startling  run 
that  ends  only  when  the  wind  subsides,  combustible  matter  is  exhausted,  or  pre¬ 
cipitation  arrives.  Normal  firebreaks,  such  as  lakes,  rivers,  or  strips  of  land 
nearly  barren  of  combustibles  may  not  terminate  spread  owing  to  branding. 

Clearly  it  is  the  wind-aiding,  not  the  mode  of  ignition,  that  is  the  key  common 
factor  in  most  fire  catastrophes. 

What  is  desired  in  analysis  of  fire  spread  (whether  in  a  single  structure 
or  on  an  urban  scale)  is  the  capacity  to  predict  with  confidence  the  rate  of 
flame  spread,  i.e.,  the  expected  position  of  the  fire  front  at  future  time,  given 
its  position  at  the  current  time.  What  is  limiting  to  meaningful  prediction  is 
quantitative  formulation  (and  then  solution)  of  the  principal  physical  processes 
controlling  spread.  With  such  Insight,  and  with  appropriate  meteorological, 
topographical,  and  fuel-loading  data,  current  computer  storage  and  speed  of 
processing  seem  adequate  to  display  fire-front  positions  at  future  time  in  the 
absence  of  countermeasures;  in  concept,  this  Lagrangian-type  calculation  is 
quite  simplistic. 

The  fastest  spread  rates  are  almost  invaviably  associated  with  wind-aiding: 
upslope  spread  usually  exceeds  downslope  spread;  spread  under  a  sustained  breeze 


usually  exceeds  spread  in  a  calm.  The  more  rapid  spread  may  be  ascribed  to 
several  factors:  hot  product  gas  blown  downwind  preheats  previously  uninvolved 
fuel  in  the  fire  path;  more-distant  transport  of  lofted  firebrands  is  likely; 
and  bent-over  plumes  may  ignite  downwind  fuel  by  contact  or  by  enhanced  radia¬ 
tive  transfer  (larger  view  factor). 

While  Interest  may  center  on  an  urban  environment  blasted  into  disarray, 
it  has  already  been  noted  that  (1)  the  debris-strewn  setting  in  the  thermo¬ 
nuclear  aftermath  has  a  far  more  continuous  distribution  of  combustible  material 
than  the  fire-code-satisfying  preblast  city,  and  (2)  the  interiors  of  (possibly 
partially  toppled)  structures  are  likely  to  be  opened.  Thus,  the  closer  proximity 
to  a  wildlands  setting  makes  it  perhaps  worthwhile  to  note  that  in  a  forest  set¬ 
ting,  flame  spread  is  usually  through  the  large-pore  matrix  of  the  understory 
(slash,  brush,  grass,  litter,  down  woody  matter),  with  an  occasional  crown  taken. 
However,  in  severe,  high-wind  conditions,  a  "wall"  of  flame  consumes  all  the 
readily  combustible  fuel  from  understory  to  overstory  (canopy)  in  one  tall  front. 
Only  very  rarely  does  flame  race  from  crown  to  crown,  either  in  the  absence  of 
an  understory  fire  or  far  in  advance  of  the  surface-level  fire  (Van  Wagner  1977)’. 
Though  the  taking  of  a  crown  is  spectacular,  aside  from  radiative  transfer  the 
event  may  not  be  that  much  more  significant  than  the  exothermicity  contributed 
by  reaction  of  a  comparable  mass  of  understory  fuel.  However,  the  particular 
point  perhaps  most  worth  noting  from  forest  fires  is  that  the  small  diameter, 
thin,  leafy  matter  (which  is  relatively  quickly  dried  out)  is  consumed  as  the 
fire  front  passes;  it  is  the  loading  of  thin  fuels  that  is  pertinent  to  the  rate 
of  front  progression.  The  thicker  fuels  are  dried  out  and  consumed  (if  ever  con¬ 
sumed  at  all)  on  a  time  span  longer  than  that  associated  with  passage  of  the  fire 
front.  There  may  be  a  flare  up  behind  the  front  as  a  certain  thickness  of  fuel 
finally  is  consumed,  but  such  flaming  usually  is  isolated  from  the  fire-front 
energetics. 

Hence,  it  may  be  useful  conceptually  to  subdivide  the  structure  of  a  wind- 
aided  fire  spread  into  (1)  a  broad  preheat  zone  in  which  previously  uninvolved 
fuels  are  desiccated  and  brought  up  to  temperatures  at  which  there  is  signifi¬ 
cant  outgassing  of  (at  least  partly  combustible)  vapor  via  heat  from  the  flaming 


zone;  (2)  the  relatively  narrow  vigorous  burning  front  itself,  a  region  of 
strong  convection  in  which  the  thermally  thin  fuels  are  consumed  rapidly;  and 
(3)  the  burn-up  zone,  in  which  thicker  fuels  are  dried  and  burned,  often  by 
heterogeneous  as  well  as  homogeneous  combustion.  Wind -aided  flame  spread 
through  a  porous,  vertically  extensive  fuel  bed  entails  (at  least)  two- 
dimensional  spread,  with  reactants  initially  in  different  phases;  furthermore, 
energetics  cannot  be  treated  too  cavalierly  because  heat  that  dries  and  gasi¬ 
fies  fuel  is  not  recovered  to  sustain  the  spread. 

The  strongly  buoyant  updraft  may  form  a  "barrier"  to  the  oncoming  flow  in 
the  two-dimensional,  line-fire  case  to  be  examined  here;  even  in  the  three- 
dimensional  case,  the  low-level  oncoming  flow  may  be  largely  diverted  around  a 
strongly  buoyant  column.  Of  course,  if  the  oncoming  wind*  is  strong  enough,  it 
should  be  able  to  blow  over  the  convective  column  (so  fire  at  the  front  is  con¬ 
fined  to  surface  fuels  only);  in  contrast  for  a  weak  enough  oncoming  flow,  the 
aspect  of  wind-aiding  would  not  be  of  consequence  (and  unless  radiation  to  down¬ 
wind  fuel  brings  that  fuel  to  pyrolysis  temperature  despite  radiative,  conduc¬ 
tive,  and  convective  heat  loss  from  the  fresh  fuel  to  the  air,  there  may  be  a 
minimum  wind  below  which  there  is  no  spread).  For  an  intermediate-strength  on 
coming  flow,  the  convective  updraft  should  remain  fairly  vertical,  and  tilted 
only  modestly  downwind.  But  there  may  be  a  tendency  for  the  configuration  to  be 
in  the  intermediate,  nearly  vertical  configuration  for  a  range  of  oncoming  flow 
speeds,  because  if  the  plume  is  blown  over,  the  ignition  by  contact  of  tne  sur¬ 
face  fuels  may  adjust  (i.e.,  enhance)  the  convective- flow  strength  to  meet  the 
oncoming-flow  strength  (Taylor  1962). 


The  wind  is  constant  neither  in  magnitude  nor  in  direction.  Hence,  utilization 
of  fire- front- propagation  insight  will  probably  entail  adopting  a  (well-justified) 
quasisteady  approximation:  the  propagation  of  flame  normal  to  the  front  depends 
only  on  the  component  of  wind  instantaneously  normal  to  the  front,  even  though 
that  wind  be  varying  in  magnitude  and  direction.  During  calms  the  fire  may  dimin¬ 
ish  in  intensity  so  flames  do  not  reach  as  high  at  the  front.  In  fact,  an  ellip¬ 
tical,  preferred-axis  fire- front  configuration  is  likely  to  arise  from  a  localized 
ignition  in  a  crosswind:  the  fire  at  the  flanks  tends  to  have  a  weak  aiding  wind 
(normal  to  the  front).  Thus,  it  is  hoped  that  results  from  the  line-fire  configu¬ 
ration  being  examined  are  applicable  locally  to  each  of  several  elliptically 
shaped  heads  that  typically  evolve  in  a  wind-aided  fire,  especially  for  multiple 
sites  of  ignition. 
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3.2  A  Conjecture  about  Wind-Aided  Flame  Spread 


The  preceding  discussion  suggests  that  the  multizone,  mul ti parameter  model 
necessary  to  describe  in  detail  the  structure  of  a  wind-aided  fire  front  may  be 
bypassed  for  some  purposes,  because  perhaps  one  parameter  is  of  foremost  importance. 

Verbally,  the  argument  may  be  stated  thus:  the  speed  of  the  crosswind 
implies  a  rate  of  entrainment  into  a  two-dimensional  plume,  if  that  oncoming  wind 
cannot  penetrate  the  plume.  But  the  entrainment  requirement  of  a  maintained  buoy¬ 
ant  column  over  a  line- type  source  is  dependent  on  the  strength  of  the  buoyant 
updraft,  according  to  now-classical  concepts  for  maintained  plumes.  In  turn,  the 
strength  of  the  updraft  is  dependent  on  the  intensity  of  the  line  source,  i.e., 
the  heat  per  length  per  time  released  (here,  by  exothermic  combustion).  And  this 
intensity  by  definition  is  simply  the  product  of  the  mass  of  thin-fuel  loading 
per  unit  area,  the  exothermic  heat  derived  per  mass  of  fuel  (after  adjustment  for 
drying  and  gasification),  and  the  speed  of  fire  propagation.  Thus,  through  inter¬ 
mediate  mechanisms,  the  flame  speed  is  a  function  of  the  crosswind  speed. 

There  are  several  refinements  to  the. above  discussion,  some  of  which  are  .now 
addressed.  For  plume-fixed  coordinates,  the  oncoming  wind  is  reduced  by  the  fire- 
propagation  speed.  Furthermore,  there  is  inflow  into  the  plume  from  the  downwind 
side  equal  to  the  fire-propagation  speed;  that  is,  there  is  a  counterflow  on  the 
downwind  side  of  the  plume  for  a  plume-fixed  observor.  Whether  entrainment  into 
the  plume  is  quite  so  predominantly  from  the  upwind  side  (since  the  oncoming  wind 
often  would  far  exceed  the  fi re-propagation  speed)  need  not  be  resolved  for  present 
purposes.*  Also,  the  classical  entrainment  concept  relating  plume-axis  updraft  to 
plume-edge  inflow  pertains  to  small  density  discrepancy  between  in-column  air  and 
surrounding  air.  For  large  density  discrepancy  between  in-plume  and  ambient  air, 
perhaps  the  relation  is  more  accurately  stated  in  terms  of  mass  flux  than  volumet¬ 
ric  flux.  Furthermore,  at  low  altitudes  (i.e.,  in  the  flaming  domain),  the  en¬ 
trainment  is  augmented  significantly  beyond  the  classical  value,  and  perhaps  another 
formulation  is  required.  However,  no  simple  refinement  of  the  classical  relation 
has  gained  widespread  acceptance;  initially  the  classical  relation  is  used  here. 


The  plume  probably  is  cellular  and  effectively  three-dimensional  at  low  levels,  so 
there  is  likely  to  be  entrainment  from  all  sides;  the  cells  merge  to  form  a  two- 
dimensional  plume  at  some  height  above  the  surface. 
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If  the  classical  entrainment  relation  is  adopted  tentatively,  it  is  self- 
consistent  to  adopt  the  Boussinesq  approximation  in  the  statement  of  the  con¬ 
servation  relations.  If  x  is  the  horizontal  coordinate,  y  is  the  vertical 
coordinate,  u  is  the  horizontal  velocity  component,  w  is  the  vertical  velocity 
component,  p  is  pressure,  p  is  density,  T  is  temperature,  and  subscript  infinity 
denotes  the  ambient  state  (for  which  the  density  is  invariant  over  the  height  of 
interest),  then  the  conservation  equations  for  mass,  momentum,  and  energy  for 
the  two-dimensional  plume  over  a  line  source  may  be  taken  in  the  following  form 
(Lee  and  Emmons  1963): 
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here,  subscripts  x  and  y  denote  partial  differentiation.  For  the  Boussinesq 
approximation. 
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For  Gaussian  profiles  for  the  dependent  variables,  the  axial  values  [W(y)  and 
R(y) ]  and  the  e-folding  distance,  or  plume  width,  b(y)  appear  as  follows: 

w(x,y)  =  W(y)  expj^-x2/b2(y)  ,  (7 
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Application  of  the  integral  method  to  (1),  (3),  and  (4),  under  (2)  and  (5)-(8), 
gives  three  ordinary  differential  equations  for  the  three  dependent  variables 
b(y)»  R(y),  W(y) .  That  is, 


w(x,y)dy  =  -  u(co.y)  -s-  (bW)  =-^-aW  , 

Jo 

unck.r  the  Taylor  entrainment  hypothesis 

-u(°°,y)  -  a w(0,y)  =  a  W(y)  , 
where  a  =  0.16  (Lee  and  Emmons  1963).  Similarly 
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For  a  dry-adiabatic  lap.-e  race  (neutrally  stable  atmosphere),  the  right-hand 
side  of  (4)  vanishes,  and 
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where  the  constant  of  integration  0  is  tt.e  full-infinite  intensity  of  the  line 
source;  hence,  by  (7)  and  (8), 
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Equations  (9),  (11),  and  \13)  admit  sol  tion  in  the  form 

b  =  (3  z  ,  W  =  y  ,  R  =  T  z”1 
where  the  constants  8,y,T  are  Identified  by 
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Thus,  the  updraft  on  the  axis  of  the  linearly  growing  plume  is  invariant  with 
height,  but  the  density  discrepancy  does  decrease  inversely  with  altitude,  the 
model  clearly  being  singular  at  the  source  z  s  0. 

Under  the  Taylor  hypothesis  that  the  crosswind  is  compatible  with  the  en¬ 
trainment  requirements,  in  view  of  (10)  and  of  the  expectation  that  the  entrain¬ 
ment  is  predominantly  from  the  upwind  side  of  the  two-dimensional  plume,  if  U 
is  the  crosswind, 

f  .  1V3 
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But  the  fire  intensity  Q  is  given  by 


Q  s  Qa  vf  ,  (1/ 

where  Q  is  the  effective  chemical  heat  per  mass  of  thin  fuel,  a  is  the  mass  of 
thin  fuel  per  area,  and  vf  is  the  fire-propagation  speed.  Hence, 


vf  ,  1 
U  "  8a2 
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or  vfoc(u3/a),  since  Q  is  relatively  invariant  for  most  fuels  of  interest.  The 
very  often-cited  general  relation  that  buoyant  strength  0  goes  as  the  cube  of 
the  crosswind  speed  U  in  the  near-critical  (i.e.,  almost  blow-over)  condition, 
as  displayed  in  (16),  may  be  traced  to  Rouse  (1947).  Taylor's  interpretation 
permits  relating  the  proportionality  factor  to  a  well-known  and  already- 
measured  empirical  constant,  the  entrainment  constant  a. 

According  to  (18)  the  flame-propagation  speed  vf  increases  as  the  exother- 
micity  Q  and  thin-fuel  loading  a  decrease,  but  increases  as  the  ambient  tempera¬ 
ture  T^  increases.  These  proportionalities  reflect  the  requirement  to  generate 
sufficient  exothermici ty  per  length  per  time  Q  to  "match"  the  crosswind  U.  The 
assumption  of  reaction  rate  sufficient  for  the  flow  rate  is  implicity  adopted; 
if  the  so-called  first  Damkohler  number  is  not  large  enough  (and  forced- 
convective  extinction  occurs),  one  of  the  bases  of  (18)  is  violated  and  the  rela¬ 
tion  does  not  hold.  Also,  if  a  strong  gust  of  wind  arises,  the  flame-propagation 
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rate  cannot  adjust  Instantaneously,  the  quasisteady  relation  (18)  is  inapplicable, 
and  the  plume  is  blown  over;  in  time,  the  propagation  speed  vf  may  increase  to 
restore  a  fairly  vertical  plume  and  (18)  holds  again.  Since  by  measurement 
(2a)  <  1,  from  (16)  one  expects  W  >  U.  Further,  if.  one  anticipates  (vf/U)  <  1, 
then 


(19) 


that  is,  for  too  strong  a  crosswind,  even  if  the  burning  rates  can  keep  pace, 
the  exothermici ty  capacity  of  the  fuel  loading  cannot,  and  the  plume  is  blown 
over. 

For  wind-aided  spread  through  a  coniferous  stand,  the  following  parametric 
values  hold  (Brian  Stocks,  private  communication): 


U  4  670  cm/s,  o  ®  0.3  g/cm2,  Q  4  1.9  x  1011  ergs/g.  (20) 

Hence;  (18)  gives  vf  ■  166  cm/s,  whereas  the  observed  value  is  110  cm/s.  In 
view  of  the  approximations  (Boussinesq  approximation,  neutrally  stable  ambient, 
line  source  of  heat,  etc.),  perhaps  a  factor-of-two  error  is  to  be  expected. 
However,  the  point  may  be  worth  making  that  if  the  entrainment  parameter 
were  ascribed  the  enhanced  value  0.2,  theory  and  experiment  would  be  in  close 
accord;  in  appendix  A,  it  is  argued  that  an  enhanced  value  of  the  entrainment 
parameter  may  be  appropriate  for  the  burning  zone  of  a  fire. 


4.  SUGGESTED  DIRECTIONS  FOR  FURTHER  RESEARCH 


4.1  Introductory  Comments 

Whereas  earlier  sections  have  presented  theoretical  models  of  types  of 
urban-scale  fire;  this  concluding  section  comments  on  experimental  testing  of 
models  to  give  them  credibility.  Since  urban-scale  fire  experiments  are 
iirmlausible,  full  validation  of  the  models  is  not  to  be  expected,  and  there 
will  remain  a  large  role  for  engineering  judgment  in  the  acceptance/ rejection 
of  predictions. 

More  explicitly,  relatively  few  of  the  dimensionless  parameters  pertinent 
to  urban-scale  free  burning  are  likely  to  take  on  the  same  value  in  feasible 
laboratory  (and  somewhat  larger-scale)  experiment;  thus,  execution  of  laboratory- 
scale  and  moderately  larger-scale  experiments  does  not  permit  one  to  predict 
definitively  by  dimensional  analysis  what  would  occur  on  the  urban  scale.  How¬ 
ever,  if  one  could  demonstrate,  by  comparison  against  a  wide  range  of  experi¬ 
mental  data,  that  a  theoretical  model  could  predict  (as  accurately  as  required 
for  the  user's  purposes)  physical  events  frpm  boundary  and  initial  conditions, 
the  model  is  given  credibility.  The  wider  the  range  with  successful  comparison, 
the  greater  the  credibility.  In  complicated,  interdisciplinary,  urban-scale 
phenomena,  the  range  of  experimental  data  normally  is  not  great  enough  to  encom¬ 
pass  the  actual  situation  of  ultimate  interest,  so  corroboration  of  the  model 
remains  incomplete. 

Of  course,  even  if  full-scale  experimentation  were  possible,  laboratory- 
scale  experimentation  usually  should  precede.  If  nothing  else,  prior  labora¬ 
tory  testing  makes  proper  instrumentation  (with  respect  to  resolution,  dynamic 
range,  etc.)  of  the  large-scale  test  more  likely.  But  laboratory-scale  experi¬ 
mentation  permits  attaining  much  data  relatively  quickly  and  relatively  inexpen¬ 
sively,  so  appreciable  parametric  variation  and  considerable  repetition  (to 
check  for  error)  is  possible.  There  is  likely  to  be  better  environmental  defi¬ 
nition,  more  extensive  and  sophisticated  diagnostic  instrumentation,  and  better 
isolation  of  constituent  components,  in  the  laboratory  than  in  some  possibly 
remote,  possibly  hostile  large-scale  environment.  Conversely,  in  large-scale 
tests,  relatively  few  data  points  are  obtained  (with  redundancy  perhaps  being 
impractical  because  of  limited  environmental  control),  and  even  these  few  points 


are  obtained  sometimes  at  large  intervals  of  time.  There  is  often  a  temptation 
to  change  too  many  parameters  from  one  test  to  the  next,  so  field  tests  may 
degenerate  Into  anecdotes  (isolated,  incompletely  monitored  events  of  uncertain 
reproducibility) . 

4.2  Brief  Description  of  Useful  Experiments 

Some  specific  experimental  tasks  that  seem  to  merit  undertaking  are  now 
very  briefly  sketched.  The  ordering  is  suggestive  of  priority. 

First,  it  has  already  been  noted  that  an  item  of  particular  uncertainty  in 
the  "nuclear-winter"  scenario  is  whether  or  not  saturation  and  precipitation 
owing  to  intrnse-fi re-induced  convection  in  a  moist  ambient  would  lead  to  early 
rainout  of  significant  amounts  of  combustion-generated  soot,  such  that  the  upper- 
tropospheric  carbon-particulate  layer  would  contain  less  mass.  Whether  soot  is 
wetted  by  the  condensed  moisture,  whether  soot  tends  to  coagulate  within  drop¬ 
lets,  whether  significant  contributions  to  the  energetics  are  connected  with 
such  phenomena,  etc.,  ought  to  be  clarified  experimentally ;  while  existing 
cloud-chamber  facilities  may  be  modified  to  undertake  the  simulation  of  the  "black- 
rain"  phenomenon,  perhaps  specially  designed  facilities  are  called  for. 

Second,  over  the  last  decade  or  so  experimental  burns  in  sections  of  the 
Canadian  National  Forest  have  been  carried  out  in  connection  with  suppression 
of  spread  of  tree  diseases;  Brian  Stocks  of  the  Great  Lakes  Forest  Research 
Station  of  the  Canadian  Department  of  the  Environment  has  participated  in  the 
planning  and  execution  of  these  tests,  which  are  often  conducted  in  heavily  fuel- 
laden  coniferous  stands  (with  total  thin-fuel  loading  of  about  0.3  g/cm2).  One 
type  of  test  typically  has  concerned  wind-aided  spread  across  a  firebreak- 
isolated  170  m  x  170  m  pine-type  stand,  with  the  entire  upwind  edge  ignited 
effectively  simultaneously;  cinematography  affords  not  only  the  rate  of  spread 
across  a  stand  of  known  loading  and  fuel  type  under  known  wind,  temperature,  and 
humidity,  but  also  provides  insight  into  smoke  formation,  stability  of  a  planar 
fire  front,  and  fire- front  structure.  Not  only  would  the  results  of  already- 
executed  tests  furnish  useful  data  with  which  to  test  theoretical  models,  but 
execution  of  further  tests  would  be  worthwhile  [for  example,  a  test  is  here  advo¬ 
cated  in  which  only  the  (say)  left-half  of  the  leading  edge  would  be  ignited,  to 


examine  edge  effects,  l.e.,  flank  effects  in  wind-aided  spread].  A  totally  dis¬ 
tinct  type  of  "test"  involves  a  large-area  burn  executed  in  the  absence  of  high 
wind;  while  such  a  burn  typically  is  "forced"  artificially  (by  continual  disper¬ 
sal  of  torches)  for  rapidity  of  execution,  it  seems  worthwhile  to  ignite  just  a 
central  area  (and  to  intervene  thenceforth  only  sparingly  to  maintain  circular 
symmetry);  the  rate  of  fire-spread  outward  against  the  induced  Inflowing  wind, 
the  plume  shape  and  height  and  soot  content,  the  nature  of  branding,  etc., 
would  be  of  interest  in  large-area  fire  modeling.  While  large  experimental 
burns  are  occasionally  carried  out  in  remote  sections  of  Australia,  the  rela¬ 
tively  close  proximity  of  the  Canadian  tests  makes  them  particularly  noteworthy. 

Third,  laboratory-scale  experiments  on  wind-aided  flame  spread  through  a  well- 
defined  matrix  of  discrete  fuel  elements  (in  a  combustion  tunnel  with  provision 
for  unconstrained  plume  dynamics)  seem  a  useful  complement  to  the  Canadian  tests 
just  noted.  These  tests  afford  opportunity  to  examine  spread  with  or  against  a 
wind  of  specified  magnitude,  to  ascertain  structure,  spread  rate,  firebreak 
dimensions,  etc.  Of  course,  without  special  provision,  radiative  transfer  cannot 
assume  the  role  in  laboratory-scale  flame  that  it  can  in  field  tests.  Since  a- 
soon-forthcoming  report  by  the  authors  is  to  be  devoted  to  the  design,  construc¬ 
tion,  instrumentation,  and  test  execution  of  such  a  free-burning-fire  facility, 
only  a  few  additional  comments  are  appended  here.  The  walls  of  such  a  tunnel 
should  be  removable  so  that  one  may  also  seek  elucidation  of  some  aspects  of 
large-area  fires  in  the  absence  of  a  strong  crosswind.  In  particular,  one  can 
check  on  the  physical  mechanism  for  augmentation  of  entrainment  (over  classical 
levels  pertinent  to  higher,  weakly  buoyant  portions  of  the  plume)  into  the  flam¬ 
ing,  low-level  portion  of  a  plume  sustained  by  burning  of  synthetic  and/or 
natural  polymers.  The  conjecture  here  is  that  the  inevitably  cellular  nature  of 
the  low-level  burning  affords  greater  fire  perimeter  and  hence  increased  entrain¬ 
ment,  up  to  the  altitude  at  which  the  cells  merge  into  one  column  (above  the 
completion-of-burning  height).  This  cellular  nature  may  be  investigated  by  com¬ 
paring  the  upflux  over  several  cellular-plastic-"pool"  configurations  of  equiva¬ 
lent  total  burning  area.  It  is  anticipated  that  one  "plastic  pool"  does  not 
entrain  as  much  ambient  air  as  an  equivalent-area  configuration  with  multiple 
pools.  [Presumably  the  constraints  of  inviscid  dynamics  (i.e.,  the  Bernoulli 


equation)  preclude  Indefinitely  large  augmentation  of  entrainment  over  classical 
levels.]  The  unconstrained  plume  above  such  a  laboratory  fire  also  permits  one 
to  measure  the  time-averaged  total  mass  upflux,  lateral  temperature  profile,  and 
lateral  updraft  profile  at  each  of  several  heights  above  the  completion-of-burning 
altitude.  These  data  can  be  used  to  ascertain  whether  entrainment  is  better 
described  by  relating  plume-to-edge  quantities  associated  with  volumetric  flux  or 
with  mass  flux;  the  distinction  is  negligible  for  the  classical  context  in  which 
in-plume  and  ambient  densities  are  of  comparable  value,  but  for  highly  exothermic 
combustion  the  distinction  is  likely  to  be  of  quantitative  consequence.  While  the 
more  appropriate  formulation  of  entrainment  may  nof  *  of  indispensable  value,  it 
seems  of  enough  quantitative  significance  to  be  well  worth  clarification. 

Fourth,  while  execution  of  an  urban-scale  burn  in  a  postblast  environment  is 
impractical,  there  have  been  over  thirteen  million-acre  burns  in  North  American 
wildlands  since  1825  (Pyne  1982).  There  also  have  been  many  smaller-scale,  but 
highly  intense,  forest  and  rural-area  fires  since  World  War  II;  these  more  recent 
events. are  often  characterized  by  sophisticated  aerial  and  ground -level  monitoring 
of  the  environment,  the  topography,  the  fuel  properties,  the  low-al t.i.tude-flow  and 
plume  behavior,  the  spread  rate,  and  the  residual  unburned  matter.  While  the  earli 
of  these  huge  uncontrolled  events  entail  less-than-ideal ,  often  incomplete,  eyewit- 
r.ess-type  accounts,  it  seems  a  worthwhile  investment  to  collect,  evaluate,  and 
summarize  this  unique  source  of  data  on  urban-scale  fire  in  heavily  loaded  dis¬ 
persals  of  fuel . 

Fifth,  a  key  tenet  of  the  proposed  interpretation  of  the  firestorm  event  is 
the  appreciable  reduction  of  entrainment  whenever  a  buoyant  column  rises  in  a 
rotating  surroundings  (such  that  the  lateral  variation  of  angular  momentum  is  non- 
de<  easing).  This  reduction  of  entrainment  permits  the  relatively  undiluted 
updraft  of  strongly  buoyant  air  to  persist  to  exceptional  altitude.  The  quali¬ 
tative  validity  of  the  concept  (specifically,  that  cyclostrophic  balance  tends 
to  suppress  radial  influx  into  a  buoyant  column)  seems  established,  but  adequate 
quantitative  work  has  yet  to  appear.  The  quantification  would  be  of  value  in  a 
variety  of  meteorological  contexts,  but  since  the  firestorm  is  such  an  exceptional 
event  in  free  burning,  very  high  priority  cannot  be  given  to  the  instrumentation 
of  a  laboratory  fi rewhirl -type  apparatus  (rotating  screen  surrounding  a  pool  fire). 
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Abstract 


We  present  here  the  results  of  a  study  of  the  buoyant  column  of  air 
which  accompanies  a  very  large  fire.  In  contrast  to  earlier  studies  (of 
smaller  fires),  we  find  it  necessary  to  compute  the  pressure  deficit  in 
the  plume  for  use  in  the  significant  modification  of  the  classical  entrain¬ 
ment  rule  appropriate  for  the  burning  region.  Results  are  given  for  various 
levels  of  the  heat  release,  atmospheric  stability,  and  lateral  size. 
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1.  INTRODUCTION 


Recent  Interest  In  atmospheric  soot  deposition  from  large  fires  has  led 
to  new,  hydrocode  estimates  (Brode  et  al.  1983;  Hasslg  &  Rosenblatt  1983; 

Small  et  al.  1983)  of  plume  heights  which  are  significantly  higher  than  those 
obtained  using  traditional  plume  theory  (Morton  et  al.  19S6).  In  view  of  this 
discrepancy,  an  attempt  is  made  to  incorporate  Into  plume  theory  the  hypothe¬ 
sized  hydrodynamic  effects  of  the  pressure  deficit  inherent  in  a  tall  buoyant 
plume  of  large  diameter. 

We  consider  here  a  fully  developed  free-burning  fire  over  a  large  area 
where  there  is  neither  a  signficant  cross  wind  nor  enough  ambient  circulation 
to  imply  the  development  of  a  fire  whirl.  We  shall  construct  a  description  of 
the  axisymmetric  convective  column  over  this  quasisteady  fire  by  using  a 
modification  of  the  well-known  integral-type  treatment  of  weakly  buoyant 
plumes  (Morton  et  al  1956;  Lee  A  Emmons  1961).  Other  extensions  of  the  tradi¬ 
tional  transversely  averaged  treatment  of  weakly  buoyant  plumes  have  been  con¬ 
sidered  (Morton  1965;  Smith  1967);  here,  however,  a  fire  of  urban  scale  is  of 
interest;  tne  radius  of  the  plume  at  low  altitudes  might  be  as  large  as  2  km, 
and  its  height  could  reach  10  km;  for  these  proportions,  the  earlier  theories 
may  not  suffice. 

In  a  completely  self-contained  theory  the  burning  rate,  and  hence  the 
vertical  heat  flux,  would  be  an  output  of  the  analysis.  In  this  study,  however, 
we  postulate  the  burning  rate  and,  on  the  hypothesis  that  the  combustion  has 
been  completed  at  an  altitude  (z  *  0)  at  which  the  entrained  air  is  a  few  times 
that  needed  for  stoichiometric  burning,  we  have  an  a  priori  inference  of  the  mass 
flux  at  that  altitude.  In  such  an  analysis,  of  course,  the  domain  lies  in 
-z^  <  z  <  z2  where  -Zj  is  the  level  at  which  the  vertical  mass  flux  is  that 
corresponding  to  the  pyrolization  rate  consistent  with  the  adopted  burning  rate; 
z2  is  the  altitude  at  which  the  vertical  velocity  becomes  zero. 

For  an  Intense  urban-scale  fire,  treatment  of  the  plume  from  near-ground 
level  to  the  height  of  stagnation  of  the  updraft  requires  provision  for  (1) 
large  lateral  density  variation,  with  its  implication  that  the  conventional 
Boussinesq  approximation  is  not  uniformly  adequate;  (2)  large  vertical  rise 


through  much  of  the  troposphere  (variation  of  ambient  density  with  altitude 
must  be  included);  (3)  a  modification  of  the  classical  entrainment  relation 
(the  large  inflow  velocities  associated  with  the  burning  zone  of  large  fires 
are  not  consistent  with  the  conventional  entrainment  law);  and  (4)  a  reten¬ 
tion  of  the  pressure  gradient  in  the  momentum  balance  for  the  burning  zone. 

These  last  two  provisions  seem  not  to  have  been  considered  previously  in  non¬ 
rotating  contexts,  but  it  does  seem  apparent  that  the  large  radial  velocities 
(0[102  ft/sec])  which  are  present  in  large  fires  could  only  be  supported  by 
pressure  deficits  of  the  order  of  several  millibars.  On  an  a  priori  basis, 
it  seemed  that  this  pressure  deficit  might  also  alleviate  somewhat  the  strong 
buoyancy  effects,  but  that  alleviation  seems  to  be  of  negligible  size  over  most 
of  the  plume  length. 

In  the  real  fire,  the  chemical  heat  release  occurs  over  a  limited  axial 
distance,  which,  in  our  analysis,  is  taken  to  occur  in  -z^  <  z  <  0.  The  dis¬ 
tance  (Zj)  is  that  required  for  the  entrainment  of  an  amount  of  air  which 
ensures  the  almost  complete  combustion  of  the  pyrolizate.  In  experiments  in 
which  the  fuel  is  inserted  in  momentum- dominated  jets,  it  is  known  that  the 
required  entrainment  is  three  times  that  required  for  stoichiometric  proportions 
(J.  E.  Broadwell,  private  communication)  but,  as  we  argue  in  a  later  section,  a 
somewhat  larger  amount  of  oxidant  may  be  needed  in  buoyancy-dominated  configu¬ 
rations. 

There  has  been  a  thorough  observational  verification  of  the  Morton-Taylor- 
Turner  hypothesis  that  air  Is  entrained  In  a  weakly  buoyant  plume  according  to 
a  simple  rule  wherein  the  far-field  radial  mass  flux  at  a  given  level  is  an 
invariant  multiple,  a,  of  a  characteristic  feature  of  the  vertical  velocity  at 
the  same  altitude.  In  the  burning  zone  where  the  flow  is  very  intense,  it 
is  highly  probably  (hut  only  anecdotal  accounts  ?.re  available)  that  the  radial 
inflow  is  more  nearly  in  keeping  with  inviscid  (infinite  Reynolds  number)  hydro¬ 
dynamics.  Accordingly,  we  will  explore  here  the  Implications  of  a  burning-zone 
entrainment  which  is  more  consistent  with  observed  fire  plumes. 

A  convenient  one-parameter  characterization  of  the  stability  of  a  dry  (fire- 
prone)  atmosphere  is  furnished  by  the  polytropic  relation.  The  polytropic  con¬ 
stant  is  related  to  the  (constant)  rate  of  thermal  decrease  with  increase  of 


altitude.  It  turns  out  that  a  simple  power  law  describes  density  decrease  with 
altitude  through  the  troposphere  to  an  accuracy  quite  satisfactory  for  present 
purposes. 

While  no  radiative  transport  Is  Included  in  our  analysis,  we  anticipate 
that,  with  minor  modifications  in  the  formulation,  it  could  be  Incorporated;  we 
also  Intend  to  extend  this  study  to  take  account  of  (a)  the  evaporation  and 
possible  re-condensation  of  moisture  initially  contained  in  the  fuel,  and  (b) 
the  condensation  of  the  water  vapor  initially  contained  in  the  ambient  air. 
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2.  ENTRAINMENT 

When  a  solid  fuel  Is  distributed  over  a  more-or-les?  horizontal  surface 
in  a  very  nonuni  form  way,  the  availability  of  pyrolizaie  also  has  a  very  non- 
uniform  distribution,  not  only  because  of  the  variations  in  fuel.  t:r  unit  area 
but  also  because  of  the  nonunlformity  in  the  distribution  of  the  radiative 
flux  (from  the  fire  above)  on  that  fuel.  Consistent  with  the  foregoi  .g  vari¬ 
ability  Is  an  Intense  nonhomogeneity  In  the  density  o*  the  gas  just  Loove  the 
surface  which  Is  fomented  by  the  rather  coarse-grained  Interpenetration  of  the 
cooler  entrained  air  with  the  pyrolizate.  That  variability  in  th®  density 
distribution  is  further  enhanced  by  the  fact  that  the  chemical  heat  rclea*. 
occurs  at  the  boundary  surfaces  (interfaces)  between  blobs  of  air  ard  blebs  of 
pyrolizate.  These  enhanced  density  contrasts  imply  a  vigorous  autocor vectior 
throughout  the  layer  -z-j  <  z  <  0  (z  *  0  lies  at  the  loosely  defined  altitude 
at  which  the  combustion  of  the  rising  gas  has  been  completed).  There  is  one 
important  aspect  of  this  particular  autoconvective  layer  which  di**er  markedly 
from  others  with  which  the  authors  are  acquainted.  In  the  convective  layer  of 
our  sun,  for  example,  the  heat  is  supplied  (radiatively)  from  below  and  the 
convection  is  necessitated  by  the  instability  of  the  configuration  which  would 
prevail  In  the  absence  of  such  convection.  In  the  solar  case,  the  governing 
parameters  are  so  extreme  that  it  Is  appropriate  to  assert  that  the  horizojta1ly 
averaged  temperature  through  that  layer  has  a  distribution  with  altitude  which 
Is  very  close  to  adiabatic.  Despite  the  fact  that,  in  our  configuration,  the 
heat  is  released  within  the  convecting  layer,  the  large  lateral  scales  and  the 
intense  density  contrasts  also  imply  a  very  vigorous  convection  with  a  horizon¬ 
tally  averaged  thermodynamic  state  whose  vertical  distribution  could  ny;  b* 
appreciably  more  stable  than  the  (neutral)  adiabatic  one*  but  which  might  be 
somewhat  less  stable.  Since  the  layer  of  interest  has  only  (at  most)  a  thick¬ 
ness  of  a  few  hundred  feet,  the/ieutral  state  Is  virtually  one  of  constant 
(horizontally  averaged)  temperature.  It  Is  known  that,  in  the  burning  region 
of  the  larger  laboratory  "floor  fires",  the  entrainment  of  air  into  the  rising 
column  of  gas  greatly  exceeds  that  which  is  postulated  by  Morton  et  al.  (1956) 

Any  greater  degree  of  stability  would  act  to  depress  the  convection. 
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and  which  Is  verified  experimentally  for  a  weakly  buoyant  nonreacting  plume. 

It  is  our  hypothesis  that  this  augmented  entrainment  Is  a  consequence  of  the 
convective  process  described  in  the  foregoing.  The  quantification  of  that 
argument  follows. 

Suppose  that  the  typical  radius  of  the  rising  column  of  gas  in  one  of 

the  cells  Is  E(z)  and  that  Its  vertical  velocity  is  w.  If  It  were  not  hemmed 

In  by  other  cells,  the  entrainment  Bu  Into  that  column  (where  u  is  radial 

velocity)  would  be  aSw  under  the  Morton  et  al.  rule.  In  this  recipe,  w  Is 

several  times  W,  the  average  upward  velocity,  but,  in  the  cellular  configuration, 

a  might  be  somewhat  depressed.  Suppose  that  the  cross-sectional  area  of  the 

fire  Is  A  a  NrmrE2  -  irb2  where  N  Is  the  number  of  cells  and  m"l(z)  is  some  mean 

of  the  ratio  of  the  updraft  area  in  one  cell  to  the  total  cell  area.  The  gross 

upflow  enhancement  ( 1 . e. ,  the  gross  entrainment  bug)  at  a  level  where  the  upflow 

In  one  cell  has  speed  w  Is  bue  =  NaBw  =  (N/m^abw,  where  subscript  e  denotes 

edge.  However,  w  must  be  at  least  of  the  order  of  mW  and  It  follows  that 

b  u.  >  H^a  b  W. 
e 

There  is  an  important  caveat  to  this  argument!  The  foregoing  microentrain¬ 
ment  can  occur  at  something  close  to  the  Morton  rate  only  if  the  pressure  deficit 
at  r  ■  b  can  support  the  external  Inflow  field.  At  the  enormous  Reynolds  numbers 
of  interest  here,  the  Bernoulli  estimate  (subscript  a  denotes  ambient) 

ue(z)  3  {2Jpa(z)  “  Pe(zj]/Pe(z)  '**  » 

would  limit  that  inflow  and,  avoiding  niceties  of  detail,  we  suggest 
ue(z)  *  TS^  a  W(z),  ^2(pa(z)  -  Pe(z))/Pa(z)J  J  , 
where  TS  means  "the  smaller  of". 

We  do  not  know  N,  of  course,  but  for  urban-scale  fires,  one  can  anticipate 
confidently  that 

ue(z)  ^TS-Jn14  a  W(z),  ^(pa(z)  -  pe(z))/pa(z)J*}s  ^2(pa(z)  -  Pe(z))/pa(z) 
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Since  the  suggestion  Is  somewhat  gratuitous,  we  shall  explore  the  conse¬ 
quences  of  (2.1)  for  several  "evaluations"  of  ue(z),  the  largest  of  which  will 
be  (2.2). 

When  a  jet  of  fuel  entrains  oxidant.  It  Is  known  that  for  cases  in  which 
kinetic  rates  are  not  at  Issue,  the  oxidation  Is  completed  at  that  downstream 
position  at  which  the  entrained  air  is  thrice  that  needed  to  provide  stoichio¬ 
metric  proportions.  For  a  fuel  whose  adiabatic  stoichiometric  flame  tempera¬ 
ture  would  be  about  2100  K  and  whose  pyrolization  temperature  is  about  600  K, 
the  requirement  that  the  temperature  change  from  the  fuel  bed  surface  to  the 
end  of  the  combustion  zone  be  small  would  imply  that  the  entrained  air  be  six 
times  the  stoichiometric  value*.  With  (1.1)  and  with  this  hypothesis  concern¬ 
ing  the  amount  of  entrained  air  in  hand,  only  an  elementary  calculation  will 
be  required  to  get  a  description  of  the  macroscopic  features  of  the  lower  part 
of  the  plume. 

Above  the  fire,  the  autoconvection  will  decay  rather  quickly  with  height 
and  the  traditional  entrainment  rule  should  be  applicable.  Actually,  of  course, 
the  autoconvection  will  begin  to  decay  somewhat  below  the  "fire-completion 
altitude"  and  the  entrainment  will  have  a  smooth  transitional  character  which 
is  approximated  only  crudely  by  the  discontinuous  picture  outlined  in  the  fore¬ 
going  account.  It  is  unlikely  that  the  prediction  of  plume  height  or  of  the 
macroscopic  feature  of  the  variation  of  the  vertical  velocity  with  altitude 
would  be  sensitive  to  this  approximation. 


This  value  is  arrived  at  by  noting  that  the  burning  of  a  simple  hydrocarbon 
with  about  fifteen  times  as  much  air  by  mass  (a  typical  stoichiometric  propor¬ 
tion)  yields  an  increase  in  temperature  of  the  fuel-air  mixture  of  roughly 
1800  K.  If  one  then  admixes  five  times  more  air,  the  increase  in  temperature 
distributed  over  the  total  mixture  is  reduced  to  300  K.  But  an  increase  of 
300  K  over  ground-level  ambient  characterizes  pyrolysis  temperatures  for  poly¬ 
meric  fuels  of  interest.  Thus  combustion  with  six  times  stoichiometric  air 
implies  a  horizontally  averaged  burning-zone  temperature  which  remains  close 
to  the  pyrolysis  temperature. 
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3.  FORMULATION 

3.1  Conservation  Equations 

The  following  Inviscld  axl symmetric  conservation  relations  are  adopted 
for  the  analysis  of  the  plume  (with  subscripts  r  and  z  denoting  partial  dif¬ 


ferentiation  with  regard  to  the  cylindrical  and  axial  coordi nates ) : 

(prw)2  +  (pru)r  =  0  ;  (3.1) 

(pruw)z  +  rpf  +  (pru*)r  -  0  ;  (3.2) 

(p  rwJ)z  +  (p  r  uw)r  +  (rp)  2  +  p  rg  -  0  ;  (3.3) 

[pwrcp(T  -  Ta)]^  +[pur  cp(T  -  T#)]^  *  P  w  r  [cp  Tfif2+  g|  -  0  .  (3.4) 


The  notation  to  be  used  is  conventional:  u  and  w  are  the  radial  and  axial  com¬ 
ponents  of  velocity;  p,  p,  and  T  are  the  gas  density,  pressure,  and  temperature; 
g  is  the  magnitude  of  the  gravitational  acceleration;  cp,  cy,  R,  and  y  are,  re¬ 
spectively,  the  heat  capacity  (taken  as  constant)  at  constant  pressure,  the 
heat  capacity  at  constant  volume,  (cp-  cy),  and  (cp/cy).  The  subscript  a  denotes 
evaluation  as  r  +  »,  i.e.,  in  the  ambient,  where  quantities  are  a  function  of 
z  only. 

Since  fire-prone  ambients  are  usually  quite  low  on  moisture,  the  dry- 
adiabatic  lapse  rate, 


is  the  limiting  case  of  interest;  exceptional,  autoconvectively  unstable  lapse 
rates  [T.  <  (-g/c_)]  are  outside  the  scope  of  interest,  which  is  concentrated 

O  P 

on  the  cases  T.  _  >  (-g/c„).  The  neutrally  stable  case  is  referred  to  often  as 

®»2  P 

the  adiabatic  atmosphere,  a  name  whose  significance  is  made  more  apparent  when 
(3.4)  is  manipulated  into  the  equivalent  form: 
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In  establishing  the  equivalence  of  (3.4}  and  (3.6).  use  was  made  of  the  hydro¬ 
static  relation  In  the  ambient. 


P 


a.z 


pag  . 


(3.7) 


and  of  the  equation  of  state, 

p  =  pRT^  pa  *  pflRTa  . 


(3.8) 


Our  analysis  will  make  use  of  the  standard  Integral  method,  and  ordinary 
differential  equations  In  z  will  emerge  from  (3.1)-(3.3)  and  (3.6).  Accord¬ 
ingly,  we  define  W(z),  Viz),  b(z),  f(z),  and  a(z)  as  follows.  For  r  <  b(z), 

w(r,z)  «  W(z),  p(r,z)  =  pa(z)  [1  -  o(z)]  , 

p(r,z)  =  pa(z)  [1  -  f(z)3  ,  -p(r,z)  ru(r,z)  =  ¥(z)  r2/b2(z)  . 

But  for  r  >  b(z) , 


w(r,z)  *  0  ,  p(r,z)  *  pfl  ,  p(r,z)  s  pa(z)  ,  -p(r,z)  ru  (r.z)  *  f(z). 


Thus,  aside  from  p(r,z)  ru(r,z),  top-hat  profiles  are  adopted. 

When  these  profiles  are  substituted  into  (3.1)-(3.4),  and  when  these  equa¬ 
tions  are  integrated  over  r  from  r  =  0  to  r  +  »,  we  obtain,  with  the  use  of  (3.7) 
and  (3.8), 


jpaWb2(l  -  f)  *  2a  pflWb  , 


pabc 


fab2  W2  pa] 


|pa  (1  -  f)  W2  b2J  -  (Pab2c)z  -  pagf  b2  -  0  , 


(3.9) 

(3.10) 


(3.11) 


(3.12) 


pa  (f  -  a)  Wb5 


Pa  ^ 


f)  WbJ 


(-f)  [ln  (p«/p.1 


In  (3.9)  and  (3.10),  the  already-discussed  classical  entrainment  relation, 

V(z)  *  a  pa(z)  b(z)  W(z)  , 
is  used,  where  a  *  0(0.1). 

Equations  (3.9),  (3.11),  and  (3.12),  with  o  a  0,  are  minor  modifications 
of  the  classical  plume  equations,  and  they  have  the  same  peculiar  structure 
near  the  top  of  the  plume.  More  explicitly,  near  z  *  (where  W  -  0),  b 
(z2  -  z)~\  W  ^  (z2  -  z}\  and  f  remains  finite  (and  nonzero)  as  (z2  -  z)  ■+  0. 
Despite  this  unphysical  local  behavior,  the  suitability  of  the  predictions  of 
that  plume  theory  is  well  established. 

When  (3.9)-(3. 12)  are  used,  the  locally  commensurate  behaviors  are 
b  *  (z2  -  z)_Js,.  W  'v  (z2-z),  f  ~  (z2  -  z),  and  a  *  (z2  -  z).  As  in  the  classi¬ 
cal  theory,  the  mass  flux  through  z2  is  not  zero,  but  that  artifact  of  the 
entrainment  law  is  no  more  reprehensible  in  this  model  than  in  the  earlier  one. 
Furthermore,  as  we  shall  see,  the  contribution  of  o  to  (3.10)  and  (3.11)  is 
negligibly  small  except  in  the  immediate  neighborhood  of  z2,  and,  without  the 
loss  of  any  meaningful  accuracy,  we  neglect  the  contribution  containing  a  in 
(3.11)  and  (3.12).  Equation  (3.10)  remains  a  recipe  for  o  from  which  we  can 
obtain  the  pressure  deficit  at  z  *  0,  and  thereby  estimate  the  entrainment  in 
the  fire  zone,  z  <  0. 


4.  ORDERS  OF  MAGNITUDE 


It  Is  useful  to  estimate  the  orders  of  magnitude  of  the  state  variables 
and  of  the  length  scales  on  which  they  vary,  both  to  justify  the  foregoing 
statements  about  the  role  (or  lack  thereof)  of  the  pressure  deficit  and  also 
to  get  an  a  priori  estimate  of  the  plume  height.  We  start  with  values  at 
z  *  0  of  dependent  variables  b  and  f,  and  note  that  the  ranges  of  Interest  are 
0  <  bQ  <  2  km  and  0  <  fQ  <  Jj,  but  the  values  of  great  interest  here  are 
bQ  *  2  km  and  fQ  *  J5.  The  notation  is  f(Q)  *  fQ,  b(0)  =  bQ. 

According  to  (3.9),  the  fractional  increase  In  pWb2  over  the  height  of 
the  phenomenon  is 

A(pWb2)  aZg 
p  W  b2  b 

When  b  is  of  order  z2,  this  fractional  change  is  small  and  (3.12)  then  implies 


■-W  • 


This  relation,  in  turn,  suggests,  if  z  is  in  km. 


1  -  f 
1  -  f- 


0[exp(z/5)]  . 


For  small  fQ,  f  vanishes  at  about  z  a  5fQ;  in  particular,  for  fQ  =  h*  f  becomes 
zero  at  z  -  3.5  km.  Thus,  z2  *  0( 5fQ  km). 

The  vertical  velocity  is  driven  by  the  buoyancy  term  in  (3.11),  possibly 
with  alleviation  through  the  pressure-gradient  term,  so  W2  can  be  estimated, 
with  the  aid  of  (4.1),  by 


W2  a  0 


The  reader  may  wish  to  recall  that  there  Is  no  Implication  that  the  updraft 
W  actually  reaches  the  value  -  fj,^*5  *  160  m/sec,  but  it  does  sug¬ 

gest  that  some  moderate  fraction  of  that  value  can  be  expected.  The  largest 
fuel  loading  which  can  be  attributed  to  an  urban  area  Is  15  g/cm2.  Pyroliza- 
tion  of  that  fuel  will  occur  at  about  600  K.  If  the  horizontally  averaged  tem¬ 
perature  is  to  change  very  little  with  height,  then  (from  Section  2)  the 
entrained  air  at  the  fire-completion  level  will  be  six  times  stoichiometric. 

The  time  scale  on  which  the  preponderance  of  fuel  was  consumed  during  the  most 
intense  phase  of  the  Hamburg,  Germany  fire  of  July  28,  1943  was  about  two  hours. 
Use  of  these  numbers  gives  a  value  for  the  updraft  W  as  follows: 


W  'v 


[15(gfue1/cma)3[90  gqas/gfuel] 
[7200  sec][0. 5x10" 3  ggas/cm33 


(4.5) 


a  value  which  is  a  very  small  fraction  of  the  160  m/s  mentioned  above.  We 
shall  return  to  this  difficulty. 

The  general  size  of  the  pressure  deficit  can  now  be  estimated  using  (3.10) 
and  (4.4);  if  b(0)  *  bQ,  YPa(0)/Pa(0)  *  a0a,  then 


a 


=  0 


bo  f0  ~] 

3  aQ2  (1  -  fQ)J 


(4.6) 


which,  for  bQ  =  2  km,  gives  approximately  10  millibars.  This  exceeds  greatly 
the  pressure  deficit  needed  in  -z1  <  z  <  0  to  permit  inflow  rates  which  have 
been  observed  [in,  for  example,  the  Hamburg  fire  (Carrier  et  al-  1983)3,  or 
which  would  be  consistent  with  the  observation  that  combustion  is  completed  in 
a  height,  z^,  of  about  a  hundred  meters. 

The  discrepancies  between  the  estimates  given  in  (4.4)  and  (4.5),  and  the 
estimates  Inferred  from  experience,  suggest  one  more  characteristic  feature  to  be 
expected  in  a  large-area  fire.  It  is  known  that  the  strongly  wind-driven  confla¬ 
grations  known  as  crown  fires  (and  other  very  energetic  fires)  have  a  structure 
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resembling  that  of  Fig.  1.  In  the  fire  depicted  in  that  figure,  the  heat  which 
is  released  in  A-B,  and  the  still  unbumed  pyrollzate  which  may  accompany  the 
burned  gases,  are  carried  almost  horizontally  into  the  almost  vertically  rising 
plume.  In  an  approximately  axially  symmetric  configuration  with  an  Intense 
radial  Inflow,  we  should  expect  that  the  outer  reaches  of  the  burning  area  would 
have  the  same  convective  aspects,  and  that  the  burning  pyrolizate  would  be  swept 
at  low  altitude  (as  in  Fig.  2)  Into  the  base  of  a  column  of  gas  whose  radius  is 
smaller  than  that  of  the  burning  area  and  within  which  the  processes  described 
In  Section  2  would  take  place.  In  order  to  formulate  the  problem  in  a  way 
which,  without  completely  arbitrary  assignment  of  bQ,  incorporates  these  con¬ 
siderations,  we  define  two  extreme  cases.  In  the  first,  bQ  =  R,  the  radius  of 
the  fire.  It  Is  readily  established  that  the  exceptionally  rapid  necking 
shown  In  Fig.  3  would  characterize  the  solution.  Furthermore,  the  pressure 
deficit  at  z  *  0  for  that  solution  would  permit  the  inflow  needed  to  supply 
the  (six  times  stoichiometric)  mass  flux  through  z  =  0  to  be  incorporated  in  a 
distance  zx  of  some  50  m.  At  the  other  end  of  the  spectrum  of  remotely  plausi¬ 
ble  values  of  bQ  Is  that  for  which,  as  Indicated  in  Fig.  4,  b2(0)  is  zerc. 

For  that  hypothesized  value,  bQ  comes  out  to  be  440  m  and  z^  is  770  m.*  Neither 


*If  the  requirement  bz(0)  =  0  is  liberalized  slightly  to  the  requirement 
[pa(l-f)b2]z  *  0  at  z  =  0,  then  (3.9)  gives  for  f  *  h* 

Wz(0)  *  4a  W(0)/bo  . 

Also,  (3.11)  gives,  without  the  pressure-gradient  term, 

2W(0)  Wz(0)  »  g. 

Combining  these  two  rleations,  one  has 

8aW2(Q)  8aW2(0)  b * 

b  * - - =  - —  . 

g  g  V 

But,  from  the  estimate  (4.5),  if  R  Is  recalled  to  be  the  fire  radius, 


hence. 


W2(0)  b0"  s  (3.5  m/sec)2  R4 ; 
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of  these  values  for  bQ  gives  a  description  commensurate  with  observation  (e.g., 
Hamburg)  and  the  range  of  possibilities  iles  well  within  400  m  <  bQ  <  2000  m. 

If  it  turns  out  that  the  overall  plume  height  and  general  geometry  are  Insensi¬ 
tive  to  variations  in  bQ  over  this  range,  a  more  definitive  choice  for  bQ  is  not 
needed. 

Finally,  before  we  carry  out  the  analysis  appropriate  to  the  foregoing,  we 
close  the  order-of-magnitude  discussion  by  noting  that,  even  if  the  estimates 
of  (4.3)  and  (4.4)  were  not  too  large,  the  size  of  the  pressure  contributions  to 
each  of  (3.11)  and  (3.12)  could  be,  at  most,  of  the  order  of  ]%  of  the  more 
important  terms  except  near  Zg,  where,  as  noted  in  earlier  discussion,  the  struc¬ 
ture  given  by  the  model  is  unrealistic  but  harmless. 


or  R  ^  2  km,  bQ  *  440  m,  and  W(0)  =  70  m/s.  Continuity  considerations  in 


Zj  <  z  <  0  suggest,  if  h  s  -Zj  , 


2  o(0)  (Yp.(O) 

Hb„2  W(O)  -  2irb„  h  -  — s— 


°  Y  P.(0>  ’ 


and,  from  (4.6), 


aba 

a  a - 9 - =  1.3  x  10“ 3  . 

3  a.2 


Hence, 


b  W(0) 

h  =  - - — j—  =  770  m. 

2  (20) a 


v’CvV'';- 

j  ■ 


5.  FINAL  FORMULATION 


The  description  of  the  state  variables  In  0  <  z  <  Z£  are  those  profiles 
given  below  (3.8).  The  differential  equations  which  govern  W,  a,  f,  and  b 
are  those  approximations  to  (3.9)-(3.12)  which  emerge  when  the  influence  of 
the  pressure  deficit  on  the  dynamics  Is  Ignored,  viz., 

[paWb*0  “  f)]z-2apaWb  , 

[p3  (1  -  f)  W2b2]z  -  pagf  b2  *  0  , 

[PafWb2]z“  -  Pa  “  f)  Wb*  [-y-]  (pa/paT)]z  * 


Pabo  = 


(a  b2W2  p. 


(5.1) 


(5.2) 


(5.3) 


(5.4) 


The  boundary  conditions  applicable  to  the  first  three  of  these  are: 
pa(0)  f ( 0)  W{0)  b2(0)  *  Q  , 

where  $,  the  total  heat  release  rate,  Is  given  by 

0  -  QAm/t  , 

where  Q  is  the  heat  of  combustion  per  unit  mass  of  fuel  (typically  19  kJ/g), 
A(»irR2)  Is  the  b^rninr,'  area,  m  is  the  mass  of  fuel  per  unit  area  (as  much  as 
15  g/cm2  In  major  cities),  and  t  is  the  time  scale  on  which  the  mass  m  is  con¬ 
sumed  (about  2  hours  in  the  great  Hamburg  fire); 


Pa(0)  [l  -  f(0)]  W{0)  b*(0)  =  , 


(5.5) 


(5.5) 


(5.7) 


•*VA  ;  V* y-  V;  j, 


1  .*  I,  .  >.l  t  -f  s 


-  » v\.\- '  -  • ;  -  ’i  .  ■  •  .  »  '  .\  *  ,**  „■  -  /. .%  •  •• '  .  >.  j  r  \  •  wr  (  r  *  ■’  * ,  '  -r  • 


where  v  is  the  mass  of  air  which  must  be  mixed  with  a  unit  mass  of'  fuel  to  pro¬ 
duce  a  mixture  six  times  as  lean  as  stoichiometric. 

The  third  boundary  condition  requires  the  specification  of  b(Q),  i.e,, 

b(O)  =  bQ.  (5.8) 

From  Section  4,  the  range  of  bQ  of  interest  is  R  ^  b0  ^  bQ0,  where  bQ  B  bQ0 
implies  bz(0)=0.  Actually,  the  smallest  credible  bQ  is  that  for  which  the 
pressure  discrepancy  a  Is  adequate*,  in  accord  with  (2.2),  to  support  the 
inflow  into  the  fire  which  will  supply  the  mass  flux  given  by  (5.7);  that  is, 
bQ  must  be  such  that 

2Ttpa(0)  [2  pa(0)o(0)/pa(0)],s  bQ  Zi  -  .  (5.9) 

Unfortunately,  Zj  is  not  known  in  an  a  priori  basis  and  the  criterion  becomes 
a  recipe  for  z, . 

In  any  case,  we  may  adopt  an  inverse  method  of  solution  and  integrate  (5.1)- 
(5.3),  subject  to  the  boundary  conditions  (5.5),  (5.7),  and  (5.8),  and  to  the 
ambient- defining  relations  (3.7),  (3.8),  and  pa  n*  p^,  where  k  is  a  given  con- 
sta.nt  (in  general,  y~l<  k  <  1).  The  method  of  solution  is  inverse  in  that  several 
trial  values  for  the  unknown  b(0)  are  to  be  investigated  for  each  case  (fixed  Q, 

Q,  v,  k).  From  (5.4)  vhe  pressure  deficit  may  be  found  a  posteriori  throughout 
the  weakly  buoyant  domain  0  <  z  <  Zg.  We  reiterate  that  (from  Section  4)  the 
range  of  b(0)  is  restricted  by  the  fact  that  too  small  values  lead  to  a  <  0  and 
toe  large  values  lead  to  implausibly  rapid  necking  in  of  the  plume  just  above 
the  completion-of-burning  plane  z  a  0.  Actually,  it  seems  worth  noting  that 
results  could  be  found  at  ?.  =  0  for  b,  f,  and  W,  and  for  their  first  derivatives, 
as  well  as  for  o,  without  continuing  the  results  to  larger  z. 


*For  b(0)  slightly  smaller  than  b00,  W2  is  large  and  negative,  so  that,  from  (3.10), 
the  dimensionless  pressure  deficit  o(0)  becomes  negative. 


VVW\ 


Pa(l-f)(Wb2)z  ■  2pab(26pao/par  .  (5.12) 

Pa(l-f)(W2b2)z  -  pa  a(b)2  -  pagf  b2.  (5.13) 


The  entrainment  speed  used  in  (5.12)  Is  taken  In  accord  with  (2.2),  except  that 
a  constant  factor  0  has  been  included  to  permit  Inclusion  of  losses  from  invis- 
cid  values  (1  >  0,  where  0  is  assigned).  It  is  convenient  to  Introduce  the 
ambient  speed  of  sound  aQ,  where 

aQ  -  C  YPa(0)/pa(0)]Js  .  (5.14) 


(5.15) 

(5.16) 


We  nondimensionalize  by  introducing 


X 


U  * 


.  lUSL  z 

<v 


(5.17) 


We  obtain,  upon  solving  for  the  derivatives, 

(U2  +  1)(X2)S  *  2  U  X  -  X2  , 
X2(U2  +  1)  Ug  *  UX2  -  X  (U2  -  1)  . 


(5.18) 

(5.19) 


In  the  phase  plane. 


where  s  may  be  found  sequentially.  Since  z  ■  0  corresponds  to  s  B  0,  and  since 
X(0)  may  be  found  from  b(0)  Und  U(0)  from  W(0)]  by  (547),  (5.20)  and  (5,21)  may 
be  solved  subject  to  the  constraints 


U  *  U, 


X0.  s 


(5.22) 


where  UQ  denotes  U(0)  and  XQ  denotes  X(0).  The  stopping  condition  stated  above 
(5.11)  is  rephrased 


X2(-sx)  U(-Sj) 


X  *U„ 
0  0 


(5.23) 


Equation  (5.20)  is  now  examined  to  characterize  solutions  in  the  burning 
zone  (Figure  5);  the  quadrant  U  >  0,  X  >  0  is  of  physical  interest.  There  is 
a  singular  curve  that  behaves  as 

U  n,  1  +  -j-  X  +  ....  as  X  •>  0  ;  (5.24) 

Ux  as  X  *■  «>  .  (5.25) 

This  singular  curve  separates  solutions  of  physical  Interest  from  other  solu¬ 
tions.  Solutions  of  physical  Interest  have  Ux  <  3/4  at  large  X;  these  solutions 
permit  U  0  for  some  value  of  positive  X.  In  particular,  these  solutions  have 

dX/dU  =  0  at  'J  s  X/2;  thus,  if  XQ  >  U0/2,  then  X  monotonically  increases  as  U 

monotonlcally  decreases  and  hence  no  necking  occurs  in  the  burning  zone;  otner- 
wise,  X  initially  decreases  as  U  decreases  for  larger  values  of  U,  though  X 
Increases  for  small  enough  U.  [Solutions  with  Ux  >  3/4  are  not  of  physical  inter¬ 
est  because  U  >  1  for  all  X;  these  mathematical  solutions  are  probably  more 
conveniently  studied  by  examining  U(X),  i.e.,  by  inverting  (5.20),  such  that  the 
locus  XU  *  U2  -  1  defines  zero  slope  for  solutions.] 
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6.  RESULTS 


Since  very  few  fires  of  the  scale  of  Interest  here  have  occurred,  and 
fewer  still  have  been  reported,  only  limited  parameter  variation  (of  Q,  m/t, 

A,  v,  k)  is  undertaken.  As  already  noted,  the  Inverse  procedure  for  solution 
adopted  Is  to  specify  the  quantities  Q  and  the  product  A(m/T),  i.e.,  the 
quantities  Q  and  Q,  and  then  to  Infer  [by  means  of  trials  with  b(0)]  what  value 
of  (m/t)  Is  compatible  with  a  particular  value  of  A.  In  this  way,  for  example, 
it  is  possible  with  systematic  investigation  to  find  the  fire  behavior,  for 
fixed  fire  area  A,  as  the  fuel  loading  per  burn  time,  iti/t,  is  increased. 

For  brevity  of  reference,  the  following  set  of  parametric  assignments  is 
termed  nominal  (alternatively,  the  set  is  termed  the  Hamburg  case,  because  it 
is  the  authors'  best  estimate  of  values  pertinent  to  that  severe  fire): 

bQ  *  500  m,  cp  *  1.0-103  m2/s2/K  ,  g  *  9.8  m/s2  ,  p0  =  1.0-I05  N/m2, 

Q  »  1.865#107  J/kg  ,  Q  -  1.0-1012  W  ,  TQ  s  300  K, 

a  =  0.1  ,  3  *  1.0  ,  y  ■  1.4  ,  k  s  0.95  ,  v  =  135. 

In  accordance  with  the  last  paragraph,  the  quantities  Q  and  bQ  are  specified, 

and  the  quantities  A  and  (m/t)  derived,  in  the  inverse  approach  adopted.  For 

completeness,  computation  reveals  A  =  12  km2,  (m/t)  s  4.46*10~3  kg/m2/s;  indeed, 
• 

Q  and  bQ  were  chosen  (by  trial  and  error)  to  achieve  these  values.  Unless 
explicitly  stated  otherwise,  in  all  computations  to  be  reported,  parameters 
take  their  nominal  values. 

First,  generic  results  are  given  for  the  initial-value  problem  for  the 
burning  zone.  Figure  5  presents  solutions  to  (5,20)  and  (5.22)  for  UQ  *  5  for 
a  series  of  values  of  XQ.  The  singular  solution  (dashed  curve)  emanating  from 
U  *  1,  X  *  0  distinguishes  physically  interesting  solutions  (solid  curves)  from 
physically  irrelevant  solutions  (dotted  curves).  For  large  enough  XQ,  the 
dimensionless  plume  radius  X  increases  monotonically  as  the  dimensionless  plume 
updraft  U  decreases.  Figures  6  and  7  present  the  values  of  dimensionless  dis¬ 
tance  s  (below  the  completlong-of-burnlng  plane)  associated  with  the  values  of 
U,X  of  Figure  5;  these  results  are  obtained  from  (5.21)  and  (5.25). 


Figures  8-10  present  results  for  the  initial-value  problem  for  z  >  0, 
i.e.,  above  the  completion-of-burning  plane,  from  (5.1)-(5.8),  for  the  nominal 
case.  As  modest  necking  in  of  the  plume  occurs  with  altitude,  the  updraft  in¬ 
creases  from  about  10  m/s  to  about  50  m/s;  then  the  updraft  decreases  with 
Increasing  altitude,  such  that  the  plume  reaches  a  height  of  almost  10  km  (cf. 
results  cited  at  the  beginning  of  Section  1).  The  fractional  in-plume  density 
discrepancy  from  ambient,  f,  vanishes  several  kilometers  below  the  zero-updraft 
altitude,  while  the  fractional  in-plume  pressure  discrepancy  from  ambient,  a, 
vanishes  just  a  few  kilometers  above  the  completion-of-burning  plane  (the  large 
negative  gradient  of  o  at  the  plume  top  is  included  purely  for  curiosity  and  is 
an  artifact  of  a  formulation  that  is  incomplete  near  the  plume  top).  The  column 
radius  and  updraft  in  the  flame  zone  are  presented  in  Figure  11;  the  flame  height 
is  found  to  be  about  90  m;  over  this  altitude  the  updraft  rises  from  nil  (at  the 
ground)  to  about  10  m/s,  and  the  column  radius  decreases  from  about  2  km  to  0.5  km. 

Figures  12  and  13  present  results  for  plume  radius  b  and  updraft  W  for  the 
lowest  100  m  of  the  buoyant  plunte  above  the  completion-of-burning  plane,  z  *  0; 
these  results  are  for  the  nominal  Hamburg  data  (except  that  bQ  is  varied). 

Since  the  results  tend  to  merge  above  100  m,  it  is  seen  that  larger  values  of  bQ 
imply  implausibly  large  necking,  whereas  smaller  values  of  bQ  imply  rather  large 
updrafts  just  above  the  burning  zone.  Figures  14  and  15  present  associated  val¬ 
ues  for  b  and  W  in  the  burning  zone,  below  z  *  0  (solid  lines);  in  addition  to 
the  solid-line  results  for  8  *  1  (Inviscid  dynamics),  results  for  8  *  0.5  (losses 
ascribed  to  friction)  are  given  by  dashed  lines.  The  narrower  plumes  at  the 
completion-of-burning  plane  are  seen  to  be  characterized  by  appreciably  greater 
flame  heights  and  greater  sensitivity  to  frictional  loss  In  the  strongly  buoy¬ 
ant  region.  If  losses  enter  in  the  burning  zone,  then  smaller  fuel  loadings 
over  larger  areas  tend  to  give  the  same  weakly  buoyant-plume  properties  as  would 
smaller  burning  areas  with  heavier  loadings  under  purely  inviscid  dynamics,  the 
total  heat  release  rate  Q  being  held  fixed.  This  property  of  the  model  also  may 
be  noted  from  Table  1;  Table  1  also  indicates  how  other  parametric  variations 
from  nominal  alter  the  flame  zone  (see  also  Figures  16  and  17).  Particular  insen¬ 
sitivity  of  the  flame-zone  characteristics  to  atmospheric  stability,  reflected  by  k, 


and  particular  sensitivity  to  the  dilution  required  for  an  Isothermal  burning 
zone,  reflected  by  v,  are  noteworthy.  For  the  parametric  variations  collected 
at  the  top  of  Table  1  and  displayed  In  Figures  16  and  17,  It  Is  reiterated  that 
all  parameters  not  called  out  have  their  nominal  values;  In  particular,  bQ  * 

500  m. 

Virtually  exclusive  attention  has  been  concentrated  thus  far  on  total  heat 
release  rate  Q  =  1012  W;  this  exceptionally  large  value  is  of  especial  interest 
because  it  characterized  the  Hamburg  fire,  and  it  seems  rarely  (if  ever)  to  have 
been  exceeded.  Table  2  presents  results  for  other  heat  fluxes;  larger  heat 
fluxes  imply  larger  burning  areas  with  heavier  fuel  loadings  (or  shorter  burn-up 
time)  for  fixed  fuel  exothermi city— indeed,  very  large  heat-release  rate  is 
incompatible  with  too  narrow  a  plume,  and  rather  moderate  heat-release  rate  is 
incompatible  with  too  wide  a  plume,  according  to  the  model. 

While  complete  parametric  investigation  of  the  existing  model,  or  refinement 
of  the  model  (e.g.,  near  the  plume  top),  remains  to  be  undertaken,  perhaps  the 
next  step  ought  to  be  experimental  investigation  of  the  validity  of  the  statements 
utilized  to  represent  burning-zone  entrainment.  While  larger-scale  tests  ulti¬ 
mately  would  seem  required,  laboratory- scale  experiments  may  suffice  to  initiate 
validation  of  the  key  conjecture:  enhanced  fire  perimeter  in  the  inevitably 
multicellular  burning  zone  of  a  fire  results  in  enhanced  low-level  induction  of 
ambient  air  [with  the  limit  on  enhancement  (over  Morton-Taylor-Turner  values) 
expressed  by  invlscld  dynamics,  i.e.,  by  the  Bernoulli  equation]. 

The  study  reported  here  was  undertaken  largely  in  the  hope  that  one  could 
estimate  the  heights  to  which  the  plumes  of  very  large  fires  might  extend.  To 
generate  these  estimates  it  has  been  necessary  to  hypothesize  that,  even  when 
the  plume  does  not  have  a  large  aspect  ratio,  the  ideas  underlying  the  Morton- 
Taylor-Turner  buoyant- column  approach  are  still  valid;  the  arguments  never  really 
did  depend  on  the  slenderness  of  the  plume  and  there  is  little  reason  to  doubt 
the  validity  of  this  approach  for  the  extraction  of  the  large-scale  character¬ 
istics  of  large  fires.  It  also  was  necessary  to  hypothesize  that  the  entrain¬ 
ment  limitation  Implied  by  the  Bernoulli  recipe  could  be  used  in  a  very  simple 


form  for  the  altitudes  within  which  the  burning  occurs,  and  to  hypothesize  that 
a  strong  convective  mixing  was  implied  therein  by  the  heterogeneity  of  the  fuel- 
consumption  distribution.  As  It  happens,  such  results  as  plume  height  and  high- 
altitude  diameter  are  not  particularly  sensitive  to  variations  on  these  hypothe¬ 
ses,  and  the  validity  of  our  results  Is  more  nearly  assured  than  Is  the  detailed 
validity  of  the  underlying  conjectures.  The  results  that  emerge  suggest  strongly 
that,  even  for  a  fire  whose  fuel  bed  has  the  scale  and  “loading"  of  a  major  city, 
the  plume  is  confined  largely  to  the  troposphere,  except  possibly  in  cases  where 
very  humid  low  level  air  would  strongly  augment  the  lifting  process  (  a  mecha¬ 
nism  which  has  not  been  incorporated  herein),  and/or  in  cases  where  a  strong, 
meteorological  low-pressure  feature  was  already  centered  near  the  city  when  the 
fire  began  (Carrier  and  Fendell  1983). 
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Case 


Q(J/s) 


A(km2) 


m/t  (kg/m2/s) 


Hamburg 

1.865-107 

12.0 

4.46*10’3 

0  «  1.2*1012  VI 

1.865*107 

15.2 

4.22-10’3 

q  *  1.492-107  J/s 

1.492-107 

20.8 

3.22*10~3 

x  *  0.85 

1. 865*10’ 

12.0 

4.48*10-* 

v  =  90 

1.865*107 

4.4 

1.20* 10- 3 

For  8  s  1: 

b0(km)  *  2.0 

1.865-107 

15.3 

3. 50* 10" 3 

0  1.8 

13.1 

4.09-10"3 

1.6 

11.2 

4.80* 10’ 3 

1.4 

9.7 

5.54-10"3 

1.2 

8.5 

6.29-10’3 

1.0 

7.8 

6.92-10'3 

0.9 

7.8 

6.88*10’3 

0.7 

8.4 

6.37-10’3 

0.5 

12.0 

4.46*10’3 

0.4 

15.3 

3. 51 *10“ 3 

0.35 

15.3 

3.51-1G"3 

0.3 

10.5 

5.11 - 10” 3 

For  8  *  0.5: 

b-(km)  =  2.0 

1.865-107 

16.6 

3.23*10“3 

0  0.9 

11.0 

4.87-10"3 

0.7 

13.3 

4.04*10"3 

0.5 

21.7 

2.47-10"3 

0.4 

31.3 

1.71*1 0~ 3 

0.35 

40.1 

1.34 *10“ 3 

0.30 

—  - 

47.0 

1.14-10"3 

Table  2.  Parametric  Variation  for  Various  Heat  Fluxes 


Note:  Asterisk  (*)  denotes  cases  for  which  Zi  <  10  meters  such  that 
b„  «  R.  A  dashed  line  denotes  cases  for  which  no  solution 
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Figure  5.  Phase-plane  solution  for  the  dimensionless  updraft  U  as  a  function 
of  dimensionless  column  radius  X  for  the  burning  zone,  where  the  singular  solu¬ 
tion  (dashed  line)  distinguishes  physically  interesting  solutions  (solid  lines) 
from  unphysical  solutions  (dotted  lines).  The  parameter  XQ  characterizes  the 
values  of  X  at  U  s  UQ  *  5. 
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Figure  6.  For  the  results  of  Figure  5,  t.ht  dimensionless  vertical  distances 
below  the  completion-of-burnlng  plane  (s  =  0)  vs.  dimensionless  column  radius  X 
where  X  =  X.  at  s  *  0. 


Figure  7.  For  the  results  of  Figure  5,  the  dimensionless  vertical  distance  s 
below  the  completlon-of-bumlng  plane  (s  s  0)  vs.  dimensionless  column  updraft  U, 
where  U  =  UQ  ■  5  at  s  *  0. 
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Figure  8.  For  the  nominal  (Hamburg)  case,  solution  for  the  buoyant  plume,  from 
(5.1)-(5.8),  as  a  function  of  altitude  z  above  the  completion-of-burning  plane 
(z  s  0).  Here  F  is  the  mass  flux,  pg W b2 ( 1  -  f);  G  is  the  momentum  flux, 
pflW2b2(l  -  f);  H  is  the  buoyancy  flux,  pflfWb2. 
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Figure  9.  For  the  nominal  case  of  Figure  8,  the  plume  radius  b  and  updraft  W 
vs.  altitude  z  above  the  completion-of-burning  plane  (z  -  0). 
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Figure  10.  For  the  nominal  case  of  Figure  8,  the  fractional  density  discrepancy 
of  in-plume  air  from  ambient,  f,  and  the  fractional  pressure  discrepancy  of  in¬ 
plume  air  from  ambient,  a,  vs.  altitude  z  above  the  completion-of-burning  plane 
(2=0). 
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Figure  12.  For  the  nominal  case,  except  that  bQ[=  b(0)]  is  varied,  the  plume 
radius  b  vs.  altitude  z  above  the  completion-of-burning  plane  (z  *  0).  Only  the 
lowest  100  m  of  the  weakly  buoyant  plume  are  depicted. 
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Figure  13.  For  the  conditions  of  Figure  12  (i.e.,  the  nominal  case  except  that 
bo  Is  varied),  the  updraft  W  vs.  altitude  z  above  the  completion-of-burning 
plane  (z  -  0).  Only  the  lowest  100  m  of  the  weakly  buoyant  plume  are  depicted. 


Figure  14.  For  the  conditions  of  Figure  12  (i.e.,  the  nominal  case  except  that 
bQ  is  varied),  the  plume  radius  b  vs.  altitude  z  below  the  completion-of- 
burning  plane  (z  =  0).  The  solid  curves  are  for  6  s  1  (inviscid  dynamics);  the 
dashed  curves  are  for  $  *  0.5  (losses  ascribed  to  friction).  The  minimum  value 
of  z  for  which  solution  is  given  is  the  ground  plane,  where  b  s  R,  and  the  fire 
area  A  3  ttR2. 
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Figure  15.  For  the  conditions  of  Figure  12  (i.e.,  the  nominal  case  except  that 
bQ  is  varied),  the  updraft  W  vs.  altitude  z  below  the  completion-of-burning 
plane  (z  *  0).  The  solid  curves  are  for  B  B  1  (inviscid  dynamics);  the  dashed 
curves  are  for  B  *  0.5  (losses  ascribed  to  friction).  The  minimum  value  of  z 
for  which  solution  is  given  is  the  ground  plane,  where  b  =  R,  and  the  fires  area 
A  *  ttR2. 
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Figure  17.  For  the  conditions  of  Figure  16  [the  nominal  case  (including  bQ  = 
500  m),  except  as  explicitly  noted],  the  updraft  W  vs.  altitude  z  below  the 
completion-of-burning  plane  (z  »  0).  The  minimum  value  of  z  for  which  solution 
is  given  is  the  ground  plane,  where  b  s  R,  and  the  fire  area  A  s  irR2.  The  case 
<  =  0.85  is  indistinguishable  from  the  nominal  (Hamburg)  case,  for  which  k  =  0. 
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ABSTRACT 


Quantitative  criteria  are  sought  for  onset  of  firestorms, 
severe  stationary  (nonpropagating)  holocausts  arising  via  merger 
of  fires  from  multiple  simultaneous  ignitions  in  a  heavily  fuel¬ 
laden  urban  environment.  Within  an  hour,  surface-level  radial  in¬ 
flow  from  all  directions  sustains  a  large-diameter  convective  col¬ 
umn  that  eventually  reaches  altitude  of  about  10  km  (e.g.,  Hamburg, 
Dresden,  Hiroshima).  As  the  firestorm  achieves  peak  intensity  (a 

couple  of  hours  after  the  ignitions),  inflow  speed:.:  are  inferred  to 

2 

attain  25-50  m/s;  typically  12  km  are  reduced  to  ashes,  before 
winds  relax  to  ambient  levels  in  six-to-nine  hours.  Here  the  fire¬ 
storm  is  interpreted  to  be  a  mesocyclone  (rotating  severe  local 
storm)  .  Even  with  exceedingly  large  heat  release  sustained  over  a 
concentrated  area,  in  the  presence  of  a  very  nearly  autoconvectively 
unstable  atmospheric  stratification,  onset  of  vigorous  swirling  on 
the  scale  of  two  hours  requires  more  than  concentration  of  circula¬ 
tion  associated  with  the  rotation  of  the  earth;  rather,  a  preexist¬ 
ing.  if  weak,  circulation  appears  necessary  for  firestorm  cyclo¬ 
genesis  . 

NOMENCLATURE 

B(z,t)  s  radius  of  preexisting  mesoscale  vortex,  m 

b0(z)  s  B(z,0),m 

-  2 
B  =  reference  acceleration  for  nonadiabatic  atmosphere,  m/s 

b(z,t)  =  e-folding  radial  distance  for  plume  variables,  m 

b^  -  "plume  radius"  at  ground  level,  m 

=  specific  heat  capacity  at  constant  pressure  for  air, 
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■  time-average  strength  of  maintained  heat  source,  W 
=  centerline  mass  flux  in  plume,  kg/s 

■  centerline  axial  momentum  flux  in  plume,  N 

2 

*  gravitational  acceleration,  m/ s 

«  centerline  energy  flux  in  plume,  W 

■  "transversely  averaged"  momentum  flux  in  plume,  J 

-  strength  of  maintained  mass  source,  kg/s 

■  pressure,  Pa 

■  swirl  associated  with  plume,  m/s 

2  2 

-  gas  constant  for  air,  m  /s  -K 

■  cylindrical  radial  coordinate,  m 

s  angular  momentum  per  mass  derived  from  earth  rotation,  m  /s 

■  temperature,  K 
m  time,  s 

■  radial  velocity  component,  m/s 

*  swirl  speed  of  preexisting  mesoscale  vortex,  n/s 
=  azimuthal  velocity  component,  m/s 

*  centerline  axial  velocity  in  plume,  m/s 

-  axial  velocity  component,  m/s 

s  axial  distance  above  ground,  m 

*  distance  from  subterranean  point  source  to  ground,  m 

■  entrainment  functional  defined  in  (16)  ,  taken  as  (39) 

=  value  of  a  in  absence  of  rotation,  0.093 

*  angular  momentum  per  mass  of  preexisting  mesoscale 

2 

vortex,  m  /s 

■  y(vR) 
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volumetric-flux  equivalent  of  E  [see  (37)3,  m  /s 

aximuthal  cylindrical-polar  coordinate 

polytropic  exponent  characterizing  ambient-atmosphere 

stability  [see  (38)] 

3 

density,  kg/m 

centerline  plume  pressure  deficit,  Pa 

strength  of  maintained  momentum  source,  N;  also,  angular 

2 

momentum  per  mass  rv,  m  /s 

parameter  appearing  in  entrainment  functional  (39),  m  /s 
component  of  angular  velocity  of  earth  locally  perpen¬ 
dicular  to  surface,  s 


SUBSCRIPTS 


a  *  ambient 

i  =  initial,  i.e.,  pertaining  to  z  =  0 


INTRODUCTION 


In  the  aftermath  of  the  atomic  bombing  of  Hiroshima  [1,2],  and 
of  the  massive  incendiary  bombing  of  Hamburg  [3-5]  and  Dresden  [6], 
particularly  virulent,  long-lived,  uncontrolled  burning  occurred 
that  had  few  if  any  recorded  precedents.  About  one-half  hour  after 
multiple  simultaneous  ignitions  (in  a  heavily  fuel-laden  urban  en¬ 
vironment),  effected  either  by  radiative-precursor  heating  and 
blast-wave  disruption  [2]  or  by  contact  with  burning  magnesium  and 
gelled  gasoline  dispersed  from  clustered  cannisters  [7],  the  fires 
merged  to  form  a  rather  uniformly  burning  area  of  many  square 
kilometers.  Whereas  the  ambient  winds  were  less  than  4  m/s,  the 
mass  fire  engendered  radially  inward  winds  at  street  level  from  all 
directions;  about  2-3  hours  after  the  initiating  bombing,  these 
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winds  reached  a  peak  intensity  of  about  20  m/s,  with  some  estimates 
by  professional  firefighters  of  50  m/s.  The  radially  inward  wind 
apparently  precluded  spread  beyond  the  initially  ignited  area, 
though  virtually  everything  combustible  within  this  region  was 
burned,  before  the  winds  subsided  to  moderate  in  speed  and  variable 
in  direction  about  six  hours  after  initiation.  A  single  huge  cen¬ 
tral  convective  column,  into  which  the  hot  product  gases  flowed, 
rose  to  about  10  km.  This  rare  nonpropagating  fire,  so  distinct 
from  more  common  ambient-wind- aided  spreads,  is  termed  a  firestorm. 
The  goal  of  the  present  investigation  is  to  delineate,  from  thermo¬ 
hydrodynamic  modeling,  quantitative  criteria  for  the  onset  of  a 
firestorm;  while  some  characterization  of  structure  and  properties 
is  inevitably  entailed,  full  detailed  description  of  the  event  at 
peak  intensity  is  not  the  prime  objective. 

In  the  interpretation  given  here,  the  term  firestorm  is  apt. 

In  a  conventional  meteorological  context,  storm  suggests  cyclonic 
wind  about  a  center  of  low  surface  pressure,  with  precipitation 
from  convectively  induced  advection  [i.e.,  from  buoyancy- caused 
ascent  and  saturation  of  warm  moist  air,  with  (1)  radial  influx 
under  continuity,  and  (2)  possible  attendant  spin-up  under  conser¬ 
vation  of  angular  momentum  associated  with  earth  rotation  or  some 
locally  enhanced  level].  Hence  a  firestorm  is  a  "heat  cyclone"  [8], 
a  mesolow  in  which  the  exothermicity  of  combustion,  as  distinguished 
from  the  condensation  of  water  vapor,  induces  free  convection.  Just 
as  firestorms  are  exceptional  fire  events,  so  mesolows  (thunder¬ 


storms  with  organized  rotation,  also  referred  to  as  tornado  cyclones 
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and  supercells)  are  uncommon  relative  to  the  total  number  of  thun¬ 
derstorms,  and  are  characterized  by  horizontal  scale  of  several 
kilometers  and  life-span  of  about  six  hours  [9],  Further,  just  as 
the  mesolow  is  characterized  by  towering  cumulonimbi  ascending 
through  the  depth  of  the  troposphere  to  the  tropopause,  so  the  fire¬ 
storm  is  characterized  by  a  convective  column  ascending  to  excep¬ 
tionally  great  height,  e.g.,  10  km  at  Hamburg. 

The  observation  at  low  altitudes  of  appreciable  radial  flux 
from  all  directions  toward  the  base  of  the  centrally  sited  convec¬ 
tive  column  corroborates,  rather  than  contradicts,  the  primarily 
rotating  nature  of  the  bulk  of  the  air  motion.  Investigation  of  the 
near-surface  inflow  layer  near  the  center  of  a  virgorously  rotating 
air  mass  over  a  fixed  flat  surface  shows  that  strong,  purely  swirl¬ 
ing  motion  is  altered  to  equally  strong,  purely  radial  influx  near 
the  ground,  though  immediately  at  the  ground  the  no-slip  constraint 
holds  [10-12]. 

The  firestorm  is  the  exceptional  event  in  that  diffusive  mecha¬ 
nisms  normally  relax  spin-up,  such  that  swirling  is  either  modest 
or  nil  [13],  Allusions  to  the  parallel  between  firestorms  and  tro¬ 
pospheric  storms  in  the  general  sense  of  strong  convection  accom¬ 
panied  by  strong  surface  winds  are  frequent.  However,  pertinence 
of  the  dynamic  characteristics  of  a  rapidly  rotating  air  mass  above 
a  relatively  fixed  flat  surface  plane  has  been  emphasized  by  Ebert 
[8],  Emmons  [14],  and  Long  [15];  Ebert  and  Emmons  suggest  that  the 
rotation  of  the  air  surrounding  the  plume  suppresses  entrainment 
such  that  the  buoyant  plume  rises  to  exceptional  altitude,  while 
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Emmons  and  Long  note  that  radial  near- surface  inflow  is  consistent 
with  rapid  higher-level  swirling.  Here,  quantitative  description 
seeking  onset  criteria  is  undertaken.  It  may  be  remarked  that  the 
well-known  propensity  for  long-range,  spotting-type  (discontinuous) 
spread  of  free-burning  fire  via  firebrands  in  the  event  of  fire- 
whirls  [16,17]  suggests  that  the  spatially  confined  character  of 
recorded  firestorms  yet  may  have  exceptions. 

FORMULATION  -  OUTLINE 

A  fully  established  convective  column  exists  over  a  maintained 
heat  source  of  known  strength  (total  enthalpy  released  per  time)  in 
an  atmosphere  of  known  stratification  (temperature  as  a  function  of 
pressure,  with  relative  humidity  taken  as  negligible  for  circum¬ 
stances  of  interest).  That  is,  a  conventional  entraining,  nonrotat¬ 
ing  plume  exists  above  a  fire  approximated  as  a  point  source  of  heat, 
without  associated  release  of  momentum  per  time  or  mass  per  time. 
Ultimately,  the  point  source  is  placed  at  a  subterranean  site,  such 
that  the  plume  width  at  ground  level  is  ascribed  observed  values; 
arbitrarily  taking  the  plume  as  entraining  (in  an  adiabatic  atmos¬ 
phere  with  ground-level  properties)  and  nonrotating,  between  subter¬ 
ranean  site  and  ground,  is  accepted  as  yielding  adequate  description 
at  ground  level.  Of  course,  the  plume  width  at  ground  level  is  an 
approximation,  as  is  the  acceptance  of  finite  mass  and  momentum  in¬ 
put  at  ground  level.  Indeed,  the  strength  of  the  heat  source  is 
also  an  estimate,  since  the  fuel  loading  (average  mass  of  combusti¬ 
ble  matter  per  area),  exothermicity  (heat  released  per  mass  of  com¬ 
bustible  matter,  with  adjustment,  if  required,  for  drying  prior  to 
burning) ,  and  fire  duration  are  known,  but  the  rate  of  heat  release 


in  time  is  not  known,  and  in  fact  the  duration  is  more  rigorously 
an  output  than  an  input.  Nevertheless,  the  procedure  avoids  speci¬ 
fying  even  more  (probably  unavailable)  detail  as  input. 

A  conventional  integral- type  plume  theory  yields  the  vertical 
structure  of  this  convective  column,  in  particular,  the  width,  up¬ 
ward  speed,  and  density  discrepancy  from  ambient  as  functions  of 
altitude  above  the  source  [18].  Typically,  the  plume  spreads  as  the 
speed  of  ascent  slows,  ultimately  to  zero,  owing  to  entrainment  of 
heavier  surrounding  air  into  the  convective  column;  even  below  the 
altitude  of  stagnation  of  ascent,  the  discrepancy  in  density  between 
plume  air  and  surrounding  air  vanishes  and  reverses,  and  above  the 
height  of  zero  discrepancy  the  validity  of  the  entrainment  concept 
is  in  doubt.  In  any  case,  this  known  solution  serves  as  an  initial 
condition  to  the  questions  of  interest  here:  as  a  function  of  heat- 
source  strength  and  dimension,  ambient  stratification,  and  ambient 
circulation,  how  long  would  spin-up  require  and  how  is  the  structure 
of  the  convective  column  altered.  If  spin-up  would  require  more 
than  a  couple  of  hours,  relaxation  processes  not  included  in  the 
formulation  are  inferred  to  enter  and  to  preclude  cyclogenesis. 

Seeking  answers  motivates  three  innovations.  The  aspect  of 
time  development  must  be  added  to  the  description  of  flow  over  a 
maintained  source  (as  distinguished  from  time  dependence  associated 
with  the  starting  plume  or  buoyant  thermal) .  The  plume  is  taken  to 
respond  in  a  quasisteady  manner  to  the  consequences  of  spin-up  of 
the  surrounding  air,  which  introduces  explicit  temporal  derivatives. 
Accordingly,  conservation  o.f  angular  momentum  enters;  conservation 
of  radial  momentum  must  be  treated  (since  existence  of  appreciable 


swirl  implies  existence  of  a  signituxant  radial  pressure  gradient), 
in  contrast  with  a  simple  statement  of  radially  invariant  pressure, 
to  which  the  conservation  of  radial  momentum  degenerates  for  a  non¬ 
rotating  column.  Finally,  the  effective  reduction  of  the  entrain- 
ment  parameter  with  spin-up  of  the  surrounding  air  must  be  quanti¬ 
fied  [19,20];  with  reduced  entrainment,  the  spreading  of  the  plume 
with  altitude  is  decreased.  There  are  two  further  considerations. 
First,  if  one  concentrates  on  the  strongly  buoyant  region  near  the 
fire,  the  lateral  discrepancy  in  density  (and  temperature)  between 
air  in  the  plume  and  surrounding  air  warrants  inclusion  [21],  In 
the  weakly  buoyant  region  further  above  the  fire,  the  density  of 
air  in  the  plume  may  be  identified  with  the  density  of  ambient  air 
surrounding  the  plume  (Boussinesq  approximation);  however,  for  a 
firestorm  the  convective  column  extends  to  such  heights  that  the 
variation  of  ambient  air  with  altitude  warrants  retention  (modified 
Boussinesq  approximation).  Second,  observations  of  intense  Austra¬ 
lian  bushfires  indicate  several  occasions,  at  peak  burning,  of 
clouds  above  the  top  of  the  convection  column  clearly  being  drawn 
down  in  the  column  [22],  However,  no  provision  is  to  be  made  for 
evolution  to  a  more  intense  stage  of  rotating  flow  incorporating 
two-cell  structure  (i.e.,  provision  for  a  downdrift,  at  least  near 
the  axis  at  the  top  of  the  column,  within  an  annular  upflux) . 

FORMULATION  -  EQUATIONS 

Consider  an  axisymmetric  treatment  in  cylindrical  polar  coordi¬ 
nates  (r,9,z),  with  corresponding  velocity  components  (u,v,w) ,  of 
a  plume  above  a  maintained  enthalpy  source  of  strength  E  (given) , 


without  associated  mass-source  strength  M  or  momentum- source 
strength  $•. 
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(1) 


The  plume  radius  b+0  at  the  source;  subscript  a  refers  to  ambient 
conditions,  i.e.,  to  conditions  outside  the  plume  where  (in  the 
absence  of  rotation)  thermodynamic  quantities  are  a  known  function 
of  altitude  z  only. 

The  flow  is  so  subsonic  that  (not  too  close  to  the  source)  the 
density,  except  for  the  gravitational  term  in  the  momentum  equation, 
is  well  approximated  by  its  ambient  value,  and  variance  of  the  tem¬ 
perature  from  ambient  is  linearly  related  to  the  variance  from  ambi¬ 
ent  of  density  (only) .  Here  the  ambient  values  vary  with  altitude 
because  of  the  tropospheric  scale  of  the  phenomena,  In  the  absence 
of  adequate  models,  no  provision  for  axially  distributed  chemical 
exothermicity  is  included. 

The  conservations  equations  for  mass,  radial  momentum,  axial  momen¬ 
tum,  and  energy  are  taken  to  be  described  adequately  by  the  following: 


(prw)z  4  (pru)r  -  0  ;  (2) 
(pruw)z  4  (pru2)r  4-  rpr  -  pv2  -  0  ;  (3) 
(prw2)z  +  (pruw)r  4  rpz  4  Prg  -  0  ;  (4) 
(pwrCpT) z  4  (purCpT)r  4  prgw  *  0  .  (5) 


Conservation  of  angular  momentum  involves  time-derivative,  advective, 
convective,  and  Coriolis  terms,  the  last  being  important  only  when 
planetary  vorticity  is  a  significant  source  of  angular  momentum  ($srv) : 
pr<t>t  +  (prw<j>)z  4  (pru<f)r  4  2ftprzu  *  0  . 


(6) 


The  following  Gaussian  description  is  adopted  for  the  dependent 


variables : 


w(r, 2, t)  =  W(z,t)  expC-r2/b2(z,t) ]  , 
pa(z)  -  p(r, z, t)  -  a(z,t)  exp[-r2/b2 (z, t) ]  , 
Pa(z)  -  p(r,z,t) 


f(z,t)  exp[-r2/b2 (x, t) ]  , 


where 


<J>(r,z,t)  =  4>la(r,z,t)  +  4>lb Cr , z , t)  +  $2(r,z,t)  ,  (10) 

where 

<|)la(r,z,t)  *  S(z,t)jl-exp[-r2/b2(z,t)]j  ,  (11) 

<}ilb(r,z,t)  *  ro(z)|l-expC-r2/E2(z,t)]J  ,  (12) 

4>2(r,z,t)  ■  q(z,  t)  [r2/b(z,t)  ]exp[-r2/b2 (z, t) ]  ,  (13) 

where  4>^a  is  a  "response"  to  the  Coriolis  term  and  is  a  "response" 

to  the  initial  condition  (in  which  a  weak  diffuse  vortex  preexists, 

the  circulation  T  of  which  is  invariant  over  time  scales  of  inter- 

o 

est,  but  the  characterizing  radius  B  of  which,  though  initially 
large  so  B  >>  b,  may  decrease  appreciably  over  the  time  span  of 
interest) .  The  contribution  $2  describes  plume-scale  redistribution 
of  angular  momentum.  The  functions  to  be  found,  W,  o,  f,  b,  B,  S, 
and  q,  are  all  anticipated  to  be  positive  near  the  source;  f  goes 
to  zero  and  changes  sign  at  a  finite  distance  above  the  source. 

Also i  the  equation  of  state  gives 

p(r,z, t)  -  p(r,z,t)  RT(r,z,t)  =>pa(z)  = pa(z)  RTa(z)  .  (14) 


Hydrostatics  suffices  for  the  ambient: 


(15) 


The  Taylor  entrainment  concept  is  introduced: 


-lim  (ru)  -  abW  ,  (16) 
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where,  in  the  absence  of  swirl,  a  -  aQ  (sO.093);  in  the  presence 
of  swirl,  a  is  a  functional  of  the  swirl-associated  variables. 

Integration  of  (1)  over  r,  from  r*0  (where  u  ■  0)  to  r+® 
[where  (16)  holds],  by  use  of  (7)  and  upon  setting  p(r,z,t)  =  p  (z) 

a. 

gives 
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Multiplying  (15)  by  r  and  subtracting  from  (4)  gives,  upon 

integration  over  all  r  and  use  of  (7) -(9)  and  setting  p+p  , 
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Equation  (5)  may  be  rewritten 
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if  P^Pa  in  all  terms,  by  (14)  and  (15),  and  since  R  =  (c  -cv)  = 
cv(  Y  -  1) , 


pawr 
Y-  1 


(20) 


Integration  over  r  and  use  of  (7),  (9),  (14),  and  the  Boussinesq 


relation  [(3T/3p)^  is  evaluated  in  the  ambient] 


T(r,z,t)-Ta(z) 
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gives 


-hr  Pa^^T 
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Integration  of  (3)  over  r  yields  zero  from  the  perfect  dif¬ 
ferential,  and  the  axial  derivative  term  yields  a  contribution  of 
0(b3)  =  0[(az)3],  since,  from  the  selfsimilar  solution  for  a  point 
source  in  an  adiabatic  nonrotating  ambient,  b  ■  oz.  But  a  -  o(l), 
and  perhaps  even  smaller  under  spin  up.  That  is,  the  plume-scale 
radial  inflow  u,  small  relative  to  other  velocity  components  in  a 
nonrotating  convective  column,  is  reduced  still  further  under  spin 
up  (for  heights  above  the  surface-inflow  layer),  according  to 
experiment  [17].  Hence,  (3)  reverts  to  the  cyclostrophic  balance 
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Execution  of  the  integrals  yields 
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In  this  context  it  may  be  helpful  to  note 
ll:expi;-x»)3*  dx  .  ta2 


1 0 
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I3(T)  = 
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I4(A) 
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tt.Z+mi,  r(l)  = 


dA2  ( l.+A)  A 


in 2,  I4(l)  -  I2  .  (27) 


Plume-associated  dependent  variables  (b,W,f,q)  are  anticipated 
to  change  during  firestorm  onset  more  slowly  in  time  than  the  scale 
of  the  preexisting  vortex  B.  Accordingly,  it  is  taken  that  <J>2  << 
where  2  (4>la  +  4>lb)  .  Substitution  of  ( 10) -  C 13)  in  (6)  is  con¬ 
strued  to  yield  three  relations: 


lb.  T  l^unb)r  =  0  • 


p^ih  +  lPru$ 
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(28) 

(29) 


[pru(4>ia  +  <j>2)]  +  [pwr(4»2  + '+i>]  -  0  . 

jprtj)^  4-  |pru(t}>i  +  ,2)]  +  ^Prw(<{>2.  +  4>2^  +  20r2puj  =  0  . 


(30) 
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Integration  of  (28)  over  r  yields 


(B2)  ^  »  -  2abW  , 


(31) 


while  (30)  gives  contributions  from  the  first  and  last  terms  only: 


2ftabW  . 


(32) 


Since  0  =  2^r/24  h" 1  at  midlatitudes  and  a  <  0(0.1),  the  growth  of 
swirl  from  the  rotation  of  the  earth  is  too  slow  to  account  for  sig- 
nificant  vorticity  on  the  scale  of  two  hours;  although  (32)  is  ten¬ 
tatively  retained,  the  decrement  of  B  in  time  from  (31)  is  of  more 


importance  to  spin-up  if  rQ  has  nonnegligible  value.  Finally,  inte¬ 


gration  of  (29)  over  r  yields,  if  p  *  p  , 
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by  transposition,  with  the  aid  of  (17) .  Only  the  product  (qb)  occurs 
The  equations  may  be  brought  into  canonical  form  by  introduction 
of  the  following  functionals: 


F  -  p  b2W,  G  =  p  b2W2,  H  =  pWb2f/g,  J  -  p  b 
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Use  of  (7)-(9),  (21),  and  (34)  in  (1)  gives  at  the  source 

F  -  —  -  0,  G  =■=  -  0,  H  s  ^--U.  E  . 
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A  solution  procedure  is  to  adopt  two  half-steps  in  a  simplis¬ 
tic  form  of  time- splitting .  The  seven  variables  b,  B,  S,  W,  cr ,  q, 
f  (or,  equivalently,  F,  G,  H,  J,  S,  a,  B)  are  taken  as  known  in  z 
at  a  given  time.  In  the  first  half  step,  B,  S,  and  o  are  updated 
via  (31),  (32),  and  (24);  in  (24),  updated  values  of  S,  B,  but  not 
of  b,  q,  are  used.  In  the  second  half  step,  (17),  (18),  (22),  and 
(33)  are  used  to  update  F,  G,  H,  J;  in  (18),  the  term  d(crb2)/dz  is 
transposed  and  evaluated  using  the  value  of  cb2  from  the  first  half 
step.  While  iteration  on  q,  b  in  o  could  be  adopted,  this  procedure 
might  be  dispensed  with.  At  time  zero,  B2  >>b2,  so  discarding  all 
effects  associated  with  swirl  to  obtain  a  classical  initial  plume 
is  taken  as  an  adequate  approximation. 

NONDIMENS I0NALIZATI0N;  AMBIENT  STRATIFICATION;  ENTRAINMENT  FUNCTIONAL 
If,  in  the  translational ly  invariant  set  presented  above,  z  =  0 
denotes  gound  level,  then  nondimens ionalization  could  be  carried 
out  in  terms  of  the  quantities  p  (0),  T_(0),  dp_(0)/dz,  and  E,  where 
it  is  recalled  that  p  (0)  » pa(0)RT„(0) ;  it  is  convenient  to  define 

3  cl  a 
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.  time2 
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The  quantity  (p  apa  )„  is  positive  in  a  dry  stable  atmosphere  (the 

a  a  Z 


case  of  prime  interest  here) ,  zero  in  a  dry  neutrally  stable  atmos 
phere  (a  special  case,  the  selfsimilar  solution  of  which  in  a  non¬ 
rotating  atmosphere  is  classical) ,  and  negative  in  a  dry  unstable 

atmosphere  (of  no  interest  here) .  A  length  scale  is  furnished  by 
e2/5g-l/5.  a  time  scale)  el/5g-3/5.  fl  gpeed  scale>  el/5B2/5.  How¬ 
ever,  the  neutral-stability,  dry-adiabatic-atmosphere  case  B  »  0 
would  require  special  treatment,  so  the  advantages  of  nondimen - 
sionalization  are  taken  to  be  outweighed  by  this  limitation. 


For  purposes  of  a  one-parameter  characterization  of  stratifi¬ 
cation,  a  so-called  polytropic  atmosphere  [23]  is  adopted  (1  *  tc  >  0 
where  k  is  given  and  piecewise  constant);  i.e.,  under  (14)  and  (15) 
adoption  of 
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In  fact,  1  > k >  1 ,  since  temperature  (aside  from  possible  local 

inversion  layers)  decreases  with  altitude.  Since  firestorm-prone 
ambients  are  typically  very  dry,  discussion  of  stability  on  a  dry 
basis  is  apropos. 

In  the  absence  of  experimental  data  beyond  the  general  guid¬ 
ance  that  entrainment  continually  decreases  as  the  air  surrounding 
the  plume  is  spun  up  above  some  threshold  value,  only  a  "schematic 
functional  is  presented  [<j>  an  assigned  constant;  T  an  assigned 
function,  with  T(0)  «  1]: 


(39) 


o(b,z,t)  -o0  Tj*(b,z,t)  u  ^(b.z.t)  -  *0J|.  -gy-*0  . 

BOUNDARY  AND  INITIAL  CONDITIONS 

At  z  =  0  for  t  >  0, 

b  -  bi>  given  ;  H  =  E,  given  .  (40) 

One  could  let  b^  vary  in  time,  and  could  require  that  only  the  time 
average  of  H  at  z  -  0  has  the  cited  value;  however,  these  quantities 
should  be  nearly  constant  over  the  interval  for  a  fluid  element  to 
rise  the  vertical  extent  of  the  plume: 

fZ 

top  n 

trise<-t^  “  [w(zi.*0]  d2l  •  (41) 

•0 

Here,  W(ztQp,t)  ■*  0,  and  W  is  anticipated  to  be  but  weakly  variant 
in  time  during  firestorm  onset.  In  the  absence  of  specific  data, 
bj^  and  H(0,t)  are  held  constant. 

The  initial  condition  at  t  =  0  is  taken  as: 

S  -  0,  so  B(z,0)  "BQ(z)  specified,  (42) 

where  the  choice  for  BQ(z)  is  guided  by  the  fact  that  BQ(z)  plays 
a  role  in  the  preexisting  weak  vortex  that  b  plays  in  the  current 
intense  plume.  Physically,  one  anticipates  (on  the  basis  of  data 
for  tropical  cyclones  and  mesocyclones)  that  the  pressure  anomaly 
(i.e.,  pressure  deficit  on  the  centerline,  from  ambient  value  at 
the  same  altitude)  decreases  with  increasing  altitude  (except  pos¬ 
sibly  at  low  altitude).  Therefore,  one  anticipates  that  may 
decrease  with  altitude  and  may  vanish  at  finite  altitude.  It  is 
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convenient  to  introduce 

rQ<z) 

v<2)  *  S^TiT  *  (43) 

SIMPLIFICATION 

As  anticipated  below  (32) ,  spin-up  associated  with  the  rota¬ 
tion  of  the  earth  is  negligible  relative  to  spin-up  associated  with 
intensification  of  a  preexisting  vortex.  (The  preexisting  circula¬ 
tion  seems  owing  to  fires  from  prior  recent  raids  at  Hamburg,  and 
to  orography  at  Hiroshima.)  Hence,  all  terms  involving  S  are 
ignored.  Furthermore,  rotation  associated  with  the  convective  col¬ 
umn  itself  is  small  during  the  firestorm-onset  interval  of  prime 
interest  here.  Accordingly,  all  terms  involving  qb  are  ignored, 
and  the  equation  for  J  is  discarded  completely.  It  may  be  confirmed 
a  posteriori  that  principal  evolution  during  spin-up  is  the  decrease 
of  the  radius  B  of  the  preexisting  mesoscale  vortex.  The  problem 
devolves  to  (35),  subject  to  (40)  at  z-0  (see  below)  and  to  (42)-(43) 
at  t  =  0.  Thus,  if  the  axial  pressure  gradient  is  retained  at  all, 

(24)  simplifies  to 

ob2  *  An 2  pQ  r2b2/B2  .  (44) 

In  fact,  it  seems  reasonable  to  expect  buoyancy  to  exceed  the 
axial  pressure  gradient  in  magnitude  during  the  spin- up  interval  of 
interest.  Accordingly,  during  this  time  interval,  b,W,f  remain 
effectively  identical  to  their  initial  (swirl- free)  values  given 
by  (35)  with  the  term  d(ob2)/dz  discarded.  For  consistency,  one 
expects  <j>  <  $  in  (39),  so  a  =  aQ. 


If  pa  is- taken  invariant  over  the  range  of  altitude  z  of  inter¬ 
est  [so  pfl  is  invariant  in  z  from  (38a) ],  if  a  -  QtQ  and  k  -  1,  and 
if  d(obz) /dz  ■  0,  then  from  the  first  three  relations  of  (35)  and 
from  (40)  one  recovers  the  known  similarity  solution  [19]: 


§-  a0(z  +  z.),  W-gl-  (zfe“ 


where 
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Although  the  procedure  incurs  spurious  finite  values  for  the  mass 
flux  F  and  the  momentum  flux  G  as  z  *  0  (investigated  below) ,  the 
above  formulae  furnish  starting  values  at  z  -  0  for  the  plume  solu¬ 
tion  at  all  times  t  during  firestorm  onset.  The  conventional  pro¬ 
cedure  of  placing  an  equivalent- strength  point  heat  source  at  the 
subterranean  site  z  -  ~z±>  with  an  adiabatic  nonrotating  atmos¬ 
phere  holding  in  -z^ < z < 0  (where  |z^|  is  small  enough  relative  to 
tropospheric  scale  height  that  p_  is  effectively  constant) ,  has 
been  taken  as  acceptable  a  procedure  as  any  facile  alternative  for 
starting  the  plume  solution  at  z*0.  Below,  the  relative  insensi¬ 
tivity  of  results  to  choice  of  z^  is  established  by  varying  W(0)  and 
f(0).  The  fact  is  that,  with  spin-up,  a  near- ground- level  inflow 
is  engendered,  and  absence  of  data  giving  temporal  and  spatial 


resolution  of  the  heat  release  at  Hamburg  precludes  establishing 
the  adequacy  of  adopting  (46) . 

Further,  because  of  the  convenience  of  a  closed-form  expres¬ 
sion,  for  the  virtually  adiabatic  atmosphere  holding  in  the  perti¬ 
nent  lower  troposphere  at  Hamburg,  for  given  and  e,  (45)  and 
(46)  are  utilized  below,  over  the  spin-up  interval,  in  (31).  How¬ 
ever,  the  variation  of  ambient  density  and  pressure  with  altitude 
for  a  not-quite-adiabatic  lower  troposphere  is  later  examined  as  a 
check;  in  this  check,  (45)  and  (46)  are  used  only  to  start  the  up¬ 
ward  integration  in  z  of  (35)  under  (38) ,  results  being  sought  with 
the  axial  pressure  gradient  term  both  retained  and  discarded. 

In  short,  if  it  is  recalled  that  4>  “  rv,  one  anticipates  that 
B  decreases  from  Bq(z)  to  b  as  time  t  increases,  such  that  v  increases 
at  radius  b  from  virtually  zero  at  t = 0  to  almost  (2/3) [Bq(z) V(z) ] /b (z) , 
or  (2/3) FQ(z) /b(z)  ,  from  (12),  (31),  and  (43) --although,  as  B  +  b, 
the  domain  of  validity  of  the  theory  probably  is  being  overextended. 

(The  model  fails  as  B  +  b  because  plume  circulation  qb  should  not  be 
ignored,  and  input  to  the  convective  column  from  a  rotating  near¬ 
surface  layer  is  not  adequately  incorporated.)  The  issue  being 
investigated  is  whether  the  stratification  as  reflected  in  k,  (more 
importantly)  the  size  of  the  intensely  convective  region  as 
reflected  in  the  parameter  b^,  and  (still  more  importantly)  the  rate 
of  enthalpy  release  as  reflected  in  the  parameter  e  are  such  that 
spin- up  from  weak  preexisting-vortex  swirl  to  about  15  m/s  is 
achieved  on  the  scale  of  two  hours,  as  reported  at  Hamburg. 

Thus,  in  the  present  interpretation,  a = aQ  during  the  early 
stages  of  spin-up,  such  that  B  is  reduced  from  B  to  b,  presumably 


in  about  two  hours.  Then,  a  is  reduced  significantly  in  appreci¬ 
ably  swirling  surroundings,  so  the  plume  ascends  to  exceptional 
altitude  [see  (39)].  While  very  preliminary  evidence  encourages 
the  concept  [19],  far  more  data  are  required.  The  increase  in 
plume  altitude  with  decrease  of  n  is  suggested  only  very  qualita¬ 
tively  below. 

RESULTS 

Nominal  values  for  problem  parameters,  believed  appropriate 
for  the  Hamburg  firestorm,  are  presented  in  Table  1.  The  basis 
for  the  more  critical  of  these  parameter  values  is  now  discussed. 

An  area  of  12  km2  with  157  kg/m2  of  combustible  material  was  re¬ 
duced  (virtually  completely)  to  ashes  in  about  6  h  [24,25];  for 
adopted  exothermicity  of  1.86*  10 4  J/g  of  combustibleis ,  E  = 

0 ( 10 1  *  ergs/s),  or,  for  plausible  ambient  reference  (ground- level) 
atmospheric  conditions,  by  (37),  e  -  0(3*1012  cm3/s) .  This  is  an 
average  value  for  e :  enhanced  values  at  earlier  time  could  be 
afforded  by  reduced  values  later.  At  the  Hamburg  airport  at  the 
time  of  the  firestorm  [24],  between  0.3  and  3.0  km,  k  -  0.96,  and 
a  higher  value  almost  surely  held  in  the  Hammerbrook  district 
itself;  above  3  km  or  so,  <  =  0.75.  A  value  of  k  somewhat 
greater  than  0,95  near  the  ground,  and  somewhat  greater  than  0.75 
at  midtropospheric  levels,  typically  is  adopted,  for  the  present, 
two-layer,  piecewise- constant  polytropic  parameterization  of  atmos- 
phereic  stability.  From  data  taken  at  the  Hamburg  airport,  the 
preexisting  vortex  is  characterized  at  ground  level  [8,  24]  by 
V(0)  =4  m/s  and  B  (0)  =8  km,  so  r  (0)  =0(3.2  *  104  m2/s)  ;  the  radius 


Neglecting  the  axial  pressure  gradient  in  the  second  of  (35) 
gives  b,  W,  f  as  functions  of  z  independently  of  the  spin-up, 
except  as  reflected  in  the  value  assigned  a,  from  (35),  (40),  and 
(42).  In  particular,  Figs.  8-10  give  the  vertical  profiles  of  the 
plume-scale  dependent  variables  for  aQ,  (aQ/10),  (aQ/100),  for 
several  states  of  atmospheric  stability  reflected  in  k.  Two-order- 
of-magnitude  decrement  in  cl  from  aQ  is  required  for  doubling  of  the 
the  plume  altitude,  so  one  infers  that  the  sensitivity  of  a  to 
<K  >  <j>0)  is  substantial.  If  a  modest  role  holds  for  the  axial  pres¬ 
sure  gradient,  the  decrement  in  time  of  the  characterizing  radius  B 
of  the  preexisting  vortex  differs  under  (35)  and  (46)  from  the  decre 
ment  in  time  calculated  under  (48) .  The  distinction  derives  from 
the  fact  that  1  in  (35),  whereas  ic  =  1  under  (48);  perhaps  even 
more  importantly,  stratification  of  the  ambient  is  included  appropri 
ately  in  (35) ,  whereas  density  is  held  constant  in  the  oversimpli¬ 
fied  (but  explicit)  relation  (48).  Results  are  given  in  Fig.  11. 

A  point  raised  earlier  just  below  (46)  concerns  the  spurious 
ground-level  mass  and  momentum  flux  incurred  by  using  (45)  and  (46) 
to  initiate  the  plume  calculation.  As  a  check,  ground- level  updraft 
W(0)  is  decreased  by  a  factor  of  four,  and  simultaneously  the  ground- 
level  normalized  density  discrepancy  f(0)  is  increased  by  a  factor 
of  four,  such  that  the  ground- level  release  of  heat  per  time  H(0) 
is  invariant  while  the  mass  per  time  F(0)  is  reduced  by  a  factor  of 
four  [see  (34a)].  It  turns  out  that  for  otherwise  nominal  assign¬ 
ment  of  values  to  problem  parameters,  such  alteration  of  W(0)  and 
f(0)  results  in  a  plume  of  height  5.1  km,  as  opposed  to  a  height  of 
4.2  km  without  such  alteration.  Thus,  imprecise  ground- level  mass 


and  momentum  fluxes  may  not  alter  overall  plume  dynamics  signifi¬ 
cantly.  While,  from  (48),  such  alteration  might  seem  to  alter  low- 
tropospheric  spin-up,  in  fact  a  compensatory  radial  pressure  gradi¬ 
ent  (not  included  here)  may  well -enter  to  restore  the  influx  to  the 
values  adopted.  This  nonhydrostatic-pressure-gradient  contribution 
for  strongly  buoyant  plumes  is  being  investigated,  and  is  to  be 
described  in  a  subsequent  publication.  Within  the  weakly  buoyant 
formulation  given  here,  calculations  confirm  that,  uniformly  in 
altitude  over  the  time  span  for  firestorm  onset,  the  axial  pressure 
gradient  makes  a  fairly  modest  contribution  to  the  balance  of  iner¬ 
tial  and  buoyant  effects  in  the  conservation  of  vertical  momentum. 

CONCLUSIONS 

Current  estimates  for  firestorm-onset  criteria  are  as  follows: 
"(1)  at  least  8  pounds  of  combustible  per  square  foot  of  fire  area, 
(2)  at  least  half  of  the  structures  in  the  area  on  fire  simultane¬ 
ously,  (3)  a  wind  of  less  than  8  miles  per  hour  at  the  time,  and 
(4)  a  minimum  burning  area  of  about  half  a  square  mile"  [2,  pp.  299- 
300].  On  the  basis  of  the  study  described  above,  the  following 
alternate  criteria  are  proposed:  (1)  a  localized  heat  release  of 
order  1019  erg/s  sustained  for  at  least  2-3  hours;  (2)  a  preexist¬ 
ing  weak  vortex  characterized  near  ground  level  by  swirl  of  4  m/s 
at  radius  8  km,  such  that  preexisting  angular  momentum  (per  unit 
mass)  near  ground  level  is  3.2 *10 4  m2/s;  (3)  absence  of  a  strong 
ambient  crosswind,  with  less  than  4  m/ s  perhaps  being  adequate 
constraint,  but  with  total  absence  being  even  more  conducive  to 
firestorm  onset;  and  (4)  a  very  nearly  dry-adiabatic  lapse  rate 
holding  for  the  lowest  few  kilometers  of  the  atmosphere.  Lower- 
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tropospheric  spin-up  to  about  20  m/s  within  2-3  h  seems  plausible 
under  such  criteria.  If  the  exothermicity  of  combustibles  is  taken 
to  be  that  of  dry  woody  matter  consumed  readily  in  forest  fires, 
which  is  1.86 *10 4  J/g  or  so,  then  the  requisite  fuel  loading  appears 
to  be  about  four  times  the  8  pounds  per  square  foot  cited  earlier, 
if  an  area  of  12  km2  is  entailed  and  the  burning  continues  at  high 
intensity  for  6  h  (as  reported  at  Hamburg) .  The  onset  of  swirling 
near  the  convective-column  edge  may  be  abrupt  in  that  it  can  rise 
from  nearly  nil  levels  to  20  m/s  or  so  within  a  half  hour.  The 
background  angular  momentum  associated  with  the  rotation  of  the 
earth  is  inadequate  for  spin-up  to  the  cited  swirl  speed  on  the 
scale  of  2  hours  or  so. 

Further  work  on  plumes  whose  base  temperatures  are  0(10 3  K) 
and  which  are  accompanied  by  significant  swirl  is  impeded  by  the 
current  absence  of  answers  to  the  questions:  (1)  is  the  entrain¬ 
ment  rate  more  properly  related  to  mass  entrained  [i.e.,  p  lim(ru)] 

“  r-H» 

per  unit  of  axial  mass  flux  [bW  lim  p],  or  to  volume  entrained 

r+0 

[lim(ru) ]  per  unit  of  axial  volume  flux  (bW) ;  (2)  by  how  much  is 
r-*«> 

the  entrainment  coefficient  reduced  by  an  increase  in  swirl. 

Answers  can  be  furnished  experimentally  only.  In  fact,  it  is  clear 
from  already  published  laboratory  experiments  on  firewhirls  [19,20] 
that  reduction  in  entrainment  with  swirl  is  highly  significant. 
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Table  1.  Nominal  Values  for  Parameters 


8. 00* 10s  cm 
1 . 23 • 10 1  cm/s2 
5.00*104  cm 
1.00-107  cm2/s2-K 
9.05*1018  erg/s 
9.80-102  cm/s2 
1.00*10®  g/cm-s2 
2.87*10®  cm2/s2-K 
3.00-102  K 


r0(°) 


9.30* 10" 2 

3.20*10®  cm2/s 

2.60-1012  cm3/s 

/ 9.90-10-1 ,  z  <  2.74* 
\8.30*10-\  z  >  2.74* 


10s  cm 
10 5  cm 


-  1.40 


Pa<°) 

dPa(0)/dz 


*  1.16 *10" 3  g/cma 

-  1.17* 10" 9  g/cm4 


Fig.  1  Profile  of  the  characteristic  radius  B  of  the  preexisting  vortex  vs. 
altitude  above  ground  level  z,  at  time  t  [ignition  at  ts0,  with  B(z,0)-Bo(z) , 
adopted  as  a  plausible  initial  profile].  The  decrement  in  B  owing  to  convec- 
tively  induced  advection  is  computed  for  a  constant-density  atmosphere;  the 
calculation  is  carried  to  the  approximate  altitude  to  which  the  ambient 
remained  very  nearly  autoconvectively  unstable  at  Hamburg.  Here  the  volu¬ 
metric  flux  of  the  heat  source  e=2.586‘1012  cm3/s.  As  In  all  succeeding 
figures,  parameter  values  are  those  of  Table  1  except  as  explicitly  noted. 

Fig.  2  Swirl  speed  v,  at  time  t 3 0  at  three  altitudes  z,  vs.  radial  distance 
from  the  axis  of  symmetry,  r,  for  the  result  of  Fig.  1.  (The  peak  swirl  speed 
occurs  at  r±1.12  B,  for  fixed  altitude  and  time,  for  the  Oseen-type  vortex 
adopted.)  Initial  swirling  is  taken  quite  small.  Here  e=2.586*1012  cm3/s. 

Fig.  3  Swirl  speed  v  at  time  t-1 . 2S  h,‘for  the  case  of  Figs.  1  and  2.  Swirl¬ 
ing  remains  modest. 

Fig.  4  Swirl  speed  v  at  time  t.=2.5  h,  for  the  case  of  Figs.  1,  2  and  3.  On¬ 
set  of  intense  swirling  occurs  within  a  span  of  one  hour. 

Fig.  5  Ground-level  values  of  the  swirl  v  and  the  preexisting-vortex  radius 
B  vs.  time  t,  for  the  radius  of  peak  swirl  (r-1.12  B).  Order-of-magnitude 
decrement  of  e  extends  spin-up  to  such  long  times  that  relaxational  mecha¬ 
nisms  presumably  preclude  firestorm  onset. 

Fig.  6  Ground-level  values  of  the  swirl  v  and  the  preexisting-vortex  radius 
B  vs.  time  t,  for  the  radius  of  peak  swirl  (r-1.12  B).  Enlarging  the  radius 
of  convection,  b^,  shortens  the  spin-up  time.  Here  e  =  2 . 586 *  1 0 1 2  cm3/s. 


Fig.  7  Ground-level  values  of  the  swirl  v  and  the  preexisting  vortex  radius 
B  vs.  time  t,  for  the  radius  of  peak  swirl  ( r^l .12  B).  The  initial  radius 
B0(z)  is  increased  2  km  at  all  altitudes  (denoted  B0  INC).  This  is  presented 
both  with  V(z)  decreased  so  rQ  is  unchanged  (denoted  rQ  NOM),  and  with  V(z) 
unaltered  so  rQ  is  increased  at  all  z  (denoted  V  NOM).  Modest  increase  of 
BQ  somewhat  lengthens  spin-up  time.  Here  e=2.586*1012  cmJ/s. 

Fig.  8  Plume  radius  b  vs.  altitude  above  ground  level  z.  Ground-level  val¬ 
ues  are  obtained  from  (45)-(46) ;  deviation  of  the  axial  pressure  gradient  from 
hydrostatic  owing  to  rotation  is  neglected.  Reduction  by  two  orders  of  mag¬ 
nitude  of  the  entrainment  parameter  nearly  doubles  plume  height.  Also  dis¬ 
played  is  the  effect  of  altering  the  ambient  stratification  [decreasing  k 
from  its  value  for  neutrally  buoyant  condition  (k=1)  increases  stability, 
with  k=y"^  holding  for  an  isothermal  layer].  When  two  values  of  k  are  given 
for  one  curve,  the  first  pertains  to  z<2.74  km  and  the  second,  to  z>2.74  km. 
Although  other  sources  give  higher  values,  according  to  [24]  the  plume  reached 
only  7  km. 

Fig.  9  Vertical  component  of  velocity  on  the  axis  of  symmetry,  W,  vs.  alti¬ 
tude  above  ground  z,  associated  with  the  result  of  Fig.  8. 

Fig.  10  Centerline  density  discrepancy  from  ambient,  normalized  by  the 
ambient  density  and  multiplied  by  the  magnitude  of  the  gravitational  accelera¬ 
tion,  f,  vs.  altitude  above  ground  z,  associated  with  the  result  of  Figs.  8 
and  9.  Validity  of  the  theory  for  altitudes  for  which  f<0  is  problematic. 
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Fig.  11  Preexisting- vortex  radius  B  vs.  altitude  above  ground  z  at  several 
times  since  ignition  t,  under  nominal  parameters  (Table  1);  also,  the  pre¬ 
existing  angular  momentum  rQ  vs.  z.  Also,  swirl  speed  v(1.12B,0,t)  vs.  t. 
Here,  a  variable-density,  two-layer,  polytropic  atmosphere  is  adopted,  whereas 
a  constant-density  (neutrally  stable)  atmosphere  is  adopted  in  Fig.  4. 
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